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Abstract 

This  thesis  expires  the  usefulness  of  a  computer  aided  control  design  software 
package  called  QDES  in  determining  the  limits  of  achievable  aircraft  performance  as  it 
relates  to  controller  capability.  Modem  aircraft,  particularly  filters,  are  being  designed 
to  be  statically  unstable  to  enhance  their  maneuverability  and  perfcnmance.  It  is  possiUe 
that  the  aircraft,  although  physically  capable  of  a  certain  level  of  performance  due  to  its 
engine/airfiratne  combination,  may  be  uncontrollable  up  to  diis  level.  This  study  sought  to 
develop  a  methodology  to  use  QDES  to  make  a  preliminary  analysis  of  an  aircraft  design 
to  determine  if  there  exists  a  controller  that  will  enable  this  design  to  meet  its  performance 
goals,  and  if  not,  to  determine  what  the  limits  of  achievable  performance  are. 
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EXPLORING  QDES  AS  A  TOOL  FOR  DETERMINING  LIMITS  OF 


ACHIEVABLE  PERFORMANCE 
IN  AIRCRAFT  DESIGN 


I.  Introduction 


1.1  Background 

Modem  aircraft,  particularly  military  fighters,  are  capable  of  some  very  impressive 
levels  of  performance.  The  engine  and  airfirame  technology  has  produced  aircraft  that 
seem  to  be  limited  only  by  the  ability  of  humans  to  fly  them.  Another  limitation  can  arise 
in  the  area  of  controller  design.  In  other  words,  a  given  aircraft  may  have  the  physical 
capability  to  pull  nine  g's  or  roll  at  720  degrees  per  second,  but  the  question  to  be  asked  is, 
can  a  controller  be  found  that  will  enable  it  to  perform  at  this  level  while  maintaining 
desired  stability  margins  and  response  characteristics  such  as  rise  time,  settling  time, 
overshoot,  etc.? 

A  computer  program  called  QDES,  written  by  Dr.  Stephen  Boyd  of  Stanford 
University,  may  help  to  answer  this  questitxi.  This  program  uses  the  theory  of 
parameterization  of  all  stabilizing  controllers  in  the  parameter  Q  along  with  the  fact  that 
many  closed  loop  system  petfomance  requirements  may  be  cast  as  convex  ccmstraints  c« 
the  closed  loop  transfer  fimctiai  QDES  uses  a  convex  problem  solver  to  search  for 
a  control  law  in  the  parameter  Q  which  satisfies  the  convex  constraints.  If  one  cannot  be 


1-1 


found.  QDES  returns  a  message  to  this  effect  This  thesis  exfdores  the  possibility  of  using 
QDES  to  answer  the  crucial  question  posed  in  the  previous  paragraph. 

7.2  Objective 

The  objectives  of  this  thesis  are  two-fold.  The  first  objective  is  to  exploe  the 
feasibility  of  using  QDES  to  determine  the  limits  of  achievable  performance  of  an  aircraft 
as  it  relates  to  controller  design.  The  second  objective  is  to  derive  a  simplified  {x-ocess 
whereby  QDES,  along  with  several  associated  MATLAB  script  files,  may  be  used  by 
stxneone  not  intimately  familiar  with  aircraft  ccmtrol  to  determine  the  limits  of  achievable 
performance  of  a  proposed  aircraft  design.  All  diat  would  be  required  are  the  state  space 
representatiems  of  the  linearized  aircraft  equations  of  motion  for  a  given  flight  ccxidition. 

1.3  Assumptions 

The  following  assumptions  were  made  caiceming  this  thesis: 

•  The  aircraft  responses  to  control  surface  (feflections  are  accurately  modeled  by  a 
given  set  of  linearized  equations  of  notiai  for  a  specified  flight  ccxidition.  Per  Etkin 
[9],  this  representaticxi  makes  the  following  assumptions: 

•  The  aircraft  is  a  rigid  bexly  whose  mass  is  unchanging  with  time. 

•  The  earth  is  considered  flat. 

•  Altitude  is  constant  and  thus  atmospheric  density  and  the  acceleration  of 
gravity  are  ccxistant. 

•  The  eanh's  rotation  is  neglected. 

•  Longitudinal  equations  are  in  the  wind  axis  and  the  lateral  are  in  the  body  axis. 
These  equation  sets  are  not  coupled. 
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•  The  aircraft  has  a  plane  of  symmetry  and  the  momenta  of  spinning  rotors  is 
neglected. 

•  Hie  trimmed  angle  of  attack  is  known  fen*  die  specified  flight  conditiem. 

•  The  aircraft  surface  deflection  and  surface  deflection  rate  limits  as  well  as  the 
acceleration  and  angle  of  attack  limits  are  known. 

•  All  transfer  functions  and  transfer  function  matrices  in  this  thesis  are  real,  raticxial, 
and  linear  time  invariant.  This  is  a  requirement  for  the  use  of  QDES. 

•  Aircraft  pitch  rate  (q)  and  normal  acceleration  (N^)  can  be  measured  perfectly  and 
are  available  for  feedback.  Thus,  sensor  noise  will  not  be  ermsidered  in  the  system 
model. 

•  The  aircraft  is  considered  to  be  flying  in  still  air  and  not  subject  to  atmospheric 
disturbances.  Thus,  process  noise  will  not  be  considered  in  the  system  model. 

•  All  feedback  is  considered  positive.  This  is  ccxisistent  with  most  of  the  post  modem 
control  theory  literature. 

1.4  Scope 

The  scope  of  this  thesis  will  be  limited  to  the  following  tasks: 

•  Set  up  QDES  <mi  the  AFIT  computers  and  ensure  its  proper  operation. 

•  Use  QDES  to  explOTc  the  limits  of  performance  for  two  different  fighter  aircraft. 
This  will  be  limited  to  determining  maximum  ntamal  acceleration  (N^)  for  the 
longitudinal  mode.  These  limits  will  be  subject  to  constraints  on  the  time  response 
of  the  closed  loop  system  as  well  as  ctmstraints  on  control  usage,  system  stability 
margins,  and  angle  of  attack  (AO A)  limit^ons. 

•  Develop  a  step-by-step  procedure  for  use  of  QDES  to  ’mplement  the  process 
described  in  the  previous  paragraph  given  die  assumptions  in  Section  (1.3).  This 
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will  include  the  writing  of  any  necessary  MATLAB  script  files  to  assist  the  non¬ 
controls  oriented  engineer  in  easily  carrying  out  this  procedure. 

•  Evaluate  the  value  and  limitations  of  using  the  procedure  as  develq)ed. 

1.5  Approach 

The  following  approach  wiU  be  taken  to  accomplish  the  tasks  outlined  in  section 

(1.4): 

•  Take  the  state  space  representation  of  the  linearized  equations  of  moticm  for  the 
A-4D  aircraft  at  Mach  0.6,  15,000  ft  and  design  an  inner  loop  pitch  stability 
augmentaticMi  system  (iiuier  loop  pitch  SAS)  ccmtroller  using  root  locus  design 
techniques  for  normal  acceleration  at  the  c.g.  (N2)  to  ncffmal  acceleration  ccxnmand 
(Nzcmd)* 

•  Use  QDES  to  design  a  controller  for  the  same  problem,  to  N2cind  ^ind  compare 
it  with  the  classical  design.  This  will  build  familiarity  with  the  QDES  design  process 
and  provide  a  cross  check  of  its  results. 

•  Use  QDES  to  explore  the  limits  of  performance  (maximum  N^)  subject  to 
constraints  as  listed  in  Section  (1.4). 

•  Take  the  state  space  representaticm  of  the  linearized  equations  of  motion  for  the 
longitudinal  mode  of  the  F- 16  (a  more  challenging  aircraft  to  control)  at  Mach  0.6, 
sea  level,  and  use  QDES  to  explore  the  limits  of  performance  as  described  for  the 
A-4D. 

•  Take  the  results  of  the  above  steps  and  write  the  step  by  step  procedure  described  in 
Section  (1.4). 
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1.6  Overview 


This  thesis  is  organized  into  eight  chapters  and  three  appendices.  Chapter  n  gives  a 
brief  explanation  of  the  parameterization  of  all  stabilizing  compensators  and  how  QDES 
uses  that  theory  to  design  closed  loop  systems.  Chapter  III  defines  and  explains  the 
measures  of  performance  to  be  used  in  this  thesis.  Chapter  IV  applies  the  Inner  Loop 
Pitch  SAS  classical  control  method  to  the  A-4D  to  N2cmd  problem  to  provide  a 
baseline  for  evaluaticms  of  QDES  designs.  Chapter  V  walks  through  the  QDES  design 
process  by  using  it  to  expire  the  limits  of  performance  for  the  A-4D.  Chapter  VI 
shows  the  results  of  the  same  process  applied  to  the  F-16.  Chapter  VII  provides  a  step- 
by-step  process  for  the  non-controls  engineer  to  apply  QDES  in  evaluating  whether  a 
proposed  design  can  meet  its  desired  performance  level  from  a  control  standpoint  It  also 
presents  the  lessons  learned  in  applying  QDES  to  the  limits  of  performance  problem. 
Chapter  VIII  provides  an  evaluaticxi  of  QDES'  usefulness  for  this  task,  and  some  caveats 
and  cautions  should  one  attempt  to  do  so. 
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//.  Development  of  Theory 


2.1  Background 

QDES  uses  the  theory  of  the  parameterization  of  all  stabilizing  compensators.  In  his 
thesis,  DeLaney  [6]  presents  a  very  rigorous  review  of  this  thetny  and  its  relation  to 
QDES.  Rather  than  present  a  duplication  of  that  effort,  a  somewhat  condensed  sunmary 
of  the  information  required  to  understand  the  processes  presented  in  diis  thesis  will  be 
provided. 

2 .2  Standard  Feedback  Diagram 

The  starting  point  for  any  discussion  of  control  theory  is  the  system  representation 
or  block  diagram.  Figure  2.1  is  the  standard  feedback  or  "P-K'*  diagnun.  The  following 
definitions  iiom  Boyd  [S]  describe  the  components  and  input/output  signals  of  the  closed 
loop  system. 

P:  Plant  The  system  to  be  controlled.  In  this  thesis,  the  aircraft. 

K(s):  Controller  or  Compensator.  The  algorithm  used  the  control  processor  to 
generate  the  actuator  commands  based  on  the  sensor  and  command  signals  received. 

w:  Exogenous  Inputs.  All  input  signals  to  the  model  except  the  actuator  input.  For 
an  aircraft  one  of  these  would  be  normal  accelo^tion  command  (N2(;ind)- 

z:  Regulated  Outputs.  Every  output  signal  from  the  model.  These  do  not  have  to 
be  measurable.  An  example  would  be  pitch  angle  (0). 
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exogenous  inputs  (w) 


actuator  iiyuts(u) 


PLANT 

(P) 


regulated  oytpui 


ts(z) 


saised  ougwits  (y) 


CONTROLLER 

(K(s)) 

Figure  2.1.  Standard  Feedback  Diagram 


y:  Sensed  Outputs.  Output  signals  duu  are  accessible  to  the  controller.  These  must 
be  measurable.  For  the  aircraft  Icmgitudinal  axis,  q  is  cmnmonly  used. 

u:  Actuator  Inputs.  The  inputs  to  the  model  that  can  be  manipulated  by  the 
controller.  For  most  aircraft  in  the  longitudinal  mode,  diis  is  elevator  deflection 
commanded  (Secmd)* 

The  state  space  representation  of  the  plant  transfer  matrix  P  is 


A 

Bw 

Bu 

P  = 

Dzw 

^zu 

.^y 

^yw 

Dyu 

(2.1) 


The  individual  matrices  within  P  correspond  to  those  used  in  the  state-variable  differential 
equations 


x  =  Ax-l-B^  w  +  ByU 

(2.2) 

z  C^x + D^^w+D^jjU 

(2.3) 

y  =  CyX  +  Dy^  W-H  DyuU 

(2.4) 

which  describe  the  state-space  representatioi  of  the  plant  P.  A  thorough  understanding  of 
this  plant  description  is  necessary  to  set  up  the  control  problem  in  QDES. 

Alternatively,  the  plant  matrix  can  be  partitioned  as 


Bzw 

Pzu 

Pyw 

PyU 

corresponding  to 


(2.5) 
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Z  =  +  P^uU 


y  =  PywW+PyuU 

where  Pxw  is  the  transfer  matrix  from  w  to  z,  etc. 
Finally,  the  P  matrix  can  be  partitioned  as 


(2.6) 

(2.7) 


P  = 


A 

B' 

C 

D 

(2.8) 


ctniesponding  to  the  state- variable  differential  equations 

x  =  Ax  +  Bu  (2.9) 

y  =  Cx  +  Du  (2.10) 

These  various  representations  of  the  plant  transfer  matrix  will  be  referenced 
throughout  this  thesis. 


2.3  Tranter  Matrix  Characteristics 

A  plant  transfer  matrix  as  described  in  equation  (2.8)  is  said  to  be  controllable  if 
(A,B)  is  stabilizable  and  (C^A)  is  detectable.  To  say  (A,B)  is  stabilizable  means  that  all 
right  half  s-plane  (RHP)  eigenvalues  of  A  (unstable  modes)  can  be  nx>ved  to  the  left  half 
s-plane  (LHP)  using  state  feedback.  These  modes  are  said  to  be  controllable. 
Characterizing  (C^A)  as  detectable  means  that  the  eigenvalues  of  A  which  are 
unobservable  (can't  be  seen  and  therefore  can't  be  moved)  are  in  the  open  LHP  (stable). 

To  control  a  system  it  must  be  stabilizable  and  detectable. 

A  transfer  function  or  transfer  function  matrix  is  an  element  of  RHoo  (€  RHoo)  if  it  is 

proper,  real,  rational,  and  has  no  poles  in  the  RHP  or  on  the  imaginaty  axis.  A  real, 
rational  transfer  function  can  be  written  as  a  ratio  of  finite  polynomials  in  s  with  real 
coefficients.  A  transfer  function  is  proper  if  the  degree  of  the  denominator  is  greater  than 
(M*  equal  to  the  degree  of  the  numerator  [7]. 
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2.4  Coprime  Factorization 


The  concept  of  cqmme  factorization  is  central  to  the  theory  of  the  parameterization 
of  all  stabilizing  compensators.  The  following  derivation  and  fcamulas  are  from  Francis 
[10]  as  presented  by  Ridgely  [15]. 

Two  matrices  D,  N  eRHoo  are  right  coprime  if  they  have  the  same  number  of 
columns  and  there  exist  matrices  (not  unique)  X,Y  €  RHoo  such  that 

iOD  -  YN  =  I  (Bezout  Identity)  (2.11) 

Two  matrices  D,  N  €  RHoo  are  left  cqmme  if  they  have  the  same  number  of  rows  and 

there  exist  matrices  (not  unique)  X,  Y  €  RHoo  such  that 

DX-bry'  =  I  (2.12) 

X,  Y,  X ,  and  Y  are  known  as  Bezout  Factors. 

Ridgely  [15]  shows  that  any  finite  dimensional,  linear  time-invariant  transfer  matrix  G(s) 
has  a  doubly  coprime  factorization 

G(s)  =  ND*'=  D*'n  (2.13) 

where 

which  is  known  as  a  Generalized  Bezout  Identity.  Note  that  the  (1,1)  and  (2,2)  blocks  of 
equation  (2.14)  are  equations  (2.1 1)  and  (2.12)  respectively.  This  representation  ensures 
that  G(s)  has  no  RHP  pole  zero  cancellations  guaranteeing  (A,B)  is  stabilizable  and  (C,A) 
is  detectable  as  described  in  Section  (2.3). 

To  find  the  eight  matrices  on  the  left  side  of  equation  (2.14),  start  again  with 
equation  (2.8)  with  (A,B)  stabilizable  and  (C,A)  detectable.  Choose  matrices  F  and  H 
such  that 

ApsA  +  BF  (2.15) 

AhsA  +  HC  (2.16) 


X  -yId  Yl_ri  O' 
N  dJln  xJ“[o  I_ 


(2.14) 


2-4 


are  stable.  F  can  be  a  full  state  feedback  gain  matrix  found  using  Linear  Quadratic 
Regulator  (LQR)  theory  and  H  an  estimator  gain  matrix  found  by  Linear  Quadratic 
Estimator  (LQE)  theory.  Also  deEne 

CpsC+DF  2.17) 

Bh  s  B+  HD  (2.18) 


With  equations  (2.15)-  (2. 18),  the  eight  matrices  from  equation  (2.14)  are  given  by 


Ap 

B  ■ 

N  = 

[Ah 

Cp 

D 

c 

D 

Ap 

B  ■ 

D  = 

[Ah 

H  ■ 

F 

I  J 

C 

I 

Ap 

Y  — 

[Ah 

J-Bh] 

Cp 

I 

F 

Ti 

J 

Ap 

-H' 

Y=: 

[Ah 

-H' 

F 

0 

F 

0 

(2.19) 

(2.20) 

(2.21) 

(2.22) 


2.5  Parameterization  of  All  Stabilizing  Compensators 

Consider  again  the  plant  transfer  matrix  partitioned  as  in  equation  (2.5).  Widi  this 
representation  for  P  used  in  the  block  diagram  in  Egure  (2.1),  the  closed  loop  transfer 
matrix  is 

Hzw  =  Pzw  +PzuK[l-PyuK]"^Pyw  (2.23) 

Francis  [10]  proves  that  K  stabilizes  P  if  and  only  if  it  stabilizes  Pyu*  Furthermore, 
through  block  diagram  manipulation,  it  can  be  shown  that  only  Pyy  "sees"  the  feedback 
matrix  K.  Thus,  Pyu  determines  the  stabilizability  of  the  closed  loop  system  while  the 
other  transfer  matrices  (Pzw  Pzu’  Py^)  determine  its  performance. 
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Let 


Pyu  (s)  =  ND"*  =  D'^N  (2.24) 

The  eight  RHoo  matrices  of  this  doubly-coiviine  factorizaticm  of  Pyu  can  be  found  by  the 

methods  described  in  Section  (2.4).  The  funily  of  all  real  rational  controllers  (K(s))  that 
stabilize  Pyu  and  therefore  P  is  parameterized  by 

K(s)  =  [Y (s)  - D(s)Q(s)][X(s)  -  N(s)Q(s)]~^  (right  coprime  factorization)  (2.25) 
=  [x(s)-Q(s)N(s)]  ^[y(s)-Q(s)D(s)]  Geft  coprime  factorization)  (2.26) 

where  Q(s)6  RHoo.  This  is  known  as  the  Youla  Parameterization  [15].  If  Q(s)  =  0,  K(s) 

is  called  the  Central  or  Nominal  controller.  If  matrices  F  and  H  were  found  using  LQR 
and  LQE  theory  as  in  Section  (2.4),  the  central  or  nominal  controller  is  known  as  an 
Observer  Based  controller  [10],  [15].  This  concept  is  key  to  setting  up  the  control 
problem  for  QDES. 


2 .6  Closed  Loop  Transfer  Matrices 


Starting  with  the  doubly-coprime  facti»rization  of  Pyy 


Pyu  =  ND~*=D"^N 


-Y 

D 


p  Yi_ri  O' 

N  xJ"L0  i_ 


(2.27) 

(2.28) 


and  using  the  formulas  for  K(s)  shown  in  equations  (2.25)  and  (2.26),  Francis  [10]  proves 
that  the  closed  loop  transfer  matrix  H2>y  can  be  written  as 


Hzw=Ti-T2QT3  (2.29) 

where 

Ti=P.w+PzuDYPyw  (2.30) 


T2=PzuD  (2.31) 

T3  =  DPy^  (2.32) 
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Fnuicis  also  shows  that  the  state  space  representatitHi  of  Tj,  T2,  and  T3  can  be  written  as 


Ap 

-BuF 

Bw 

Ti  = 

0 

Ah 

Bw  ■^■HDy^ 

+  D2U 

i^ZW 

T2  = 

Ap 

Bu] 

F  DzuJ 

T3  = 

Ah 

+  HDy^'j 

.^y 

- 1 

Q 

(2.33) 


(2.34) 

(2.35) 


where  the  various  matrices  are  as  defined  in  equations  (2.1)  and  (2.15)  -  (2.18).  Notice 
that  Tj,  T2,  and  T3  can  be  formed  having  mily  the  plant  infonnation  as  represented  in 
equation  (2.1)  and  its  coprime  factorization  as  described  in  section  (2.4). 


2.7  QDES  Theory  of  Operation 

QDES  seeks  to  design  the  parameter  (^s)  such  that  the  closed  loop  system  H2W 
from  equation  (2.29)  meets  a  set  of  closed  loop  constraints  or  requirements  which  are 
convex  in  Q.  Boyd  [3]  defmes  this  concept  as  follows:  "if  Hzw  is  a  closed  loop  response 
which  meets  our  perfonnance  requirements,  and  Hj^  is  any  other  closed  loop  response 
which  meets  our  requirements,  then  the  closed  loop  response  (H2w'*'H2w)/2  meets  our 
requirements."  If  the  above  condition  holds,  the  requirements  are  convex. 

QDES  is  a  C  compiler  which  translates  control  constraints  such  as  step  resptmse 
characteristics  or  transfer  function  values  into  constraints  on  Q  and  then  uses  a  numerical 
convex  problem  solver  to  find  the  values  of  Q  enabling  to  meet  these  constraints  [3]. 
If  QDES  cannot  find  the  required  values  of  Q,  this  indicates  that  no  controller  of  that 
order  can  be  designed  to  meet  them.  If  QDES  does  find  the  required  values  of  Q,  the 
controller  can  be  found  via  equations  (2.25)  and  (2.26).  Because  QDES  uses  numerical 
techniques,  it  requires  the  continuous  plant  model  to  be  discretized.  Furthermore,  QDES 
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uses  the  impulse  responses  of  the  matrices  Tj.  T2,  and  T3.  Use  of  these  in^nilse  response 
matrices  allows  the  problem  to  be  finite  dimensional  and  die  resulting  controller  (usually  of 
very  high  order)  to  be  implemented  with  feedback  around  a  multi-input,  multi-output 
(MIMO)  finite  impulse  response  (FIR)  filter  as  described  in  [3].  A  thorough  discussion  of 
the  actual  mechanics  of  using  QDES  follows  in  Ouster  VII. 

2.8  Summary 

This  chapter  discussed  the  control  theory  behind  the  operation  of  QDES.  Starting 
with  the  standard  feedback  diagram,  it  defined  the  input  and  output  vector  signals  and  the 
various  representations  of  the  transfer  matrices  that  relate  them  to  each  other.  With  these 
matrices,  coprime  factorization  was  introduced  and  was  used  to  demonstrate  the  concept 
of  the  parameterization  of  all  stabilizing  compensators.  Given  this  parameterization,  the 
closed  loop  system  Hj-vv  was  represented  in  terms  of  the  stable  closed  loop  maps  Tj,  T2, 
T3  and  the  parameter  Q.  The  concept  of  ccwivex  constraints  on  H^vv  was  then  introduced 
and  the  method  whereby  QDES  uses  a  numoical  convex  problem  solver  to  arrive  at  the 
required  values  of  Q  was  briefly  discussed. 
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III.  Measures  of  Performance 


3.1  Introduction 

As  the  puipose  of  this  thesis  is  to  explme  the  usefulness  of  QDES  in  detmnining  the 
limits  of  achievable  aircraft  performance,  it  is  necessary  to  qtecify  and  define  the  measures 
of  performance  to  be  used.  For  this  thesis,  these  measures  will  be  categmized  under 
aircraft  performance,  time  response  characteristics,  and  stability  margins. 

1.2  Aircrcft  Performance 

Aircraft  performance  can  be  characterised  by  various  measures  such  as  speed,  range, 
payload  capacity,  ordnance  delivery  accuracy,  maneuverability,  etc.  From  a  controllability 
standpoint,  maneuverability  is  the  measure  most  directly  related  to  the  aircraft  control 
system.  In  the  longitudinal  axis,  to  which  this  thesis  is  limited,  maximum  normal 
acceleration  (N2)  and  pitch  rate  (q)  are  the  two  most  closely  tied  to  the  aircraft's  control 
system.  Maximum  attainable  angle  of  attack  (AOA)  is  another  longitudinal  measure  of 
performance,  but  it  is  more  a  function  of  aerodynamic  design  rather  than  the  control 
system  capability,  unless  the  control  of  the  aircraft  at  high  AOA  is  the  area  to  be 
examined.  This  challenging  aircraft  control  problem  is  beycmd  the  scope  of  this  thesis. 
The  change  in  forward  velocir/  ;s  also  an  important  measure  of  performance,  but  it  is 
more  a  function  of  aerodynamic  design,  wing  loading,  and  thrust  availability  than  control 
system  capability. 

This  thesis  uses  N2  as  the  aircraft  measure  of  performance.  This  is  the  most  logical 
in  that  all  aircraft  are  structurally  g-limited.  Furthermore,  normal  g  capability  relates 
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directly  to  two  maneuverability  measures  inqxxtant  to  aii'to-air  combtt  and  air-to-ground 
weapons  delivery.  The  first  is  turn  radius  vt^ch  is  a  fiinction  of  true  airqieed  and  N2  . 
The  lower  the  true  airspeed  and  the  hi^ier  the  normal  g,  the  smaller  the  turn  radius  and 
the  tighter  the  turn.  The  smallest  turn  radius  occurs  at  the  well  kiK>wn  comer  velocity,  the 
lowest  airspeed  where  the  aircraft  can  |ni11  the  maximum  g  allowed  by  structural  limits 
without  exceeding  AOA  limitations.  Directly  related  to  die  turn  radius  is  the  nose  tracking 
rate,  how  many  degrees  per  second  the  nose  of  the  aircraft  moves  across  the  horizon. 
Again,  the  maximum  nose  tracking  rate  occurs  when  the  aircraft  is  pulling  maximum  g  at 
the  lowest  possible  speed. 

While  maximum  N2  relates  to  the  magnitude  of  the  aircraft's  perfmmance,  one  also 
needs  to  examine  the  quality  of  that  level  of  performance. 

3.3  Time  Response  Characteristics 

One  method  of  evaluating  the  quality  of  a  closed  loop  system’s  response,  and  thus 
the  control  law  used,  is  measurement  of  its  unit  step  response  characteristics.  QDES 
allows  easy  constraint  or  optimization  of  these  characteristics.  Furthermore,  simulation  of 
the  closed  loop  system  response  to  a  unit  step  input  provides  quick  and  easy  evaluation  of 
them. 

The  following  characteristics  and  their  definititms  wiU  be  used  in  this  thesis.  In  all 

cases,  the  command  value  is  normally  1.0,  the  unit  step  input. 

Overshoot  (Mp):  The  maximum  amount  by  which  the  step  response  exceeds  the 

command  value,  usually  expressed  as  a  percentage  of  the  command  value. 

Undershoot  (U^):  The  maximum  negative  value  (assuming  a  positive  step 
command)  attained  by  the  step  response,  usually  expressed  as  a  percentage  of  the 
command  value. 
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Rise  Tune  (Tf):  The  earliest  time  after  which  the  step  response  exceeds  0.8  of  the 
command  value  [S]. 

Settling  Time  (Tsi):  The  earliest  time  after  which  the  step  response  is  always  a 
certain  percentage  of  the  command  value  [S].  Note  that  this  is  sli^tly  different  from  the 
normally  accepted  definition  which  uses  the  final,  or  steady  state  value  rather  than  the 
command  value.  The  reason  for  this  is  that  in  using  QDES,  one  seeks  to  ccmstrain  T^ 
before  knowing  what  the  final  value  will  be.  Therefore,  it  is  assumed  to  be  the  command 
value  and  the  constraints  are  specified  acc<»dingly.  This  will  become  clear  when  the 
actual  mechanics  of  using  QDES  are  discussed  in  Chapter  V. 

Command  Error  (ec„(])‘  difference  between  the  command  value  and  the  value 
of  the  step  response  at  the  end  of  the  time  period  of  interest.  It  is  expressed  as  a 
percentage  of  the  command  value.  This  measure  is  chosen  over  the  m(»e  common  steady 
state  error  (t^)  because  in  this  thesis,  the  aircraft  response  will  be  examined  over  a  short 
time  period.  The  steady  state  error  is  defined  from  the  final  value  of  the  response  once  all 
the  transients  have  died  out.  In  the  case  of  the  Phugoid  oscillation,  this  could  be  hundreds 
of  seconds  after  the  commanded  input.  Thus,  this  measure  of  performance  would  have 
little  meaning  for  the  responses  to  be  examined. 

Figure  3. 1  illustrates  a  typical  unit  step  response. 
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It  is  not  difficult  to  appreciate  how  these  unit  step  resptmse  characteristics  relate 
directly  to  the  quality  of  the  aircraft's  response.  For  instance,  when  a  pilot  desires  to  pull 

five  g's,  he  or  she  wants  the  aircraft  to  respond  hurly  quickly  (Tf),  not  overshoot  the  mark 
excessively  (Mp),  get  to  the  commanded  value  quickly  (Tgi),  and  stay  there  (e^nid).  He  or 

she  also  doesn't  want  the  aircraft  to  initially  go  in  the  wrong  direction  (11$),  or  if  this  is 
unavoidable,  the  pilot  would  at  least  like  this  nonminimum  phase  response  to  be 
mininiized. 

3.4  Stability  Margins 

In  designing  a  control  system  for  an  aircraft  (or  any  plant),  the  control  engineer  is 
forced  to  work  with  a  mathematical  model  which  is  an  approximation  of  the  actual  system. 
The  model  ccmtains  approximation  errors  and  perhaps  unnxxleled  modes  or  dynamics. 
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The  aocumulation  of  these  errors  is  known  as  the  plant  uncoiainty.  The  control  engineer 
seeks  to  design  a  control  system  that  allows  the  closed  loop  system  to  retain  its  stability 
and  perfcnmance  characteristics  in  the  face  of  a  certain  level  of  uncertainty.  Measures  oi 
the  system's  resistance  to  being  destabilized  by  this  uncertainty  are  known  as  the  stal^ty 
margins.  Qassically,  these  are  the  gain  and  phase  margins,  that  amount  of  gain  or  i^rase 
change  that  results  in  a  change  in  the  encirclenient  of  the  critical  pc^t  (-1  for  negative 
feedback)  in  the  well  known  Nyquist  Stability  Criteria. 

A  slightly  different  representation  of  stal^ty  margins,  Imown  as  vector  margins  can 
also  be  calculated  from  the  polar  plot  of  the  c^n  loop  or  loop  transfer  function  L(s)=PK. 
The  point  of  closest  approach  of  the  polar  plot  to  the  critical  point  is  known  as  the  M 
circle  radius.  This  is  illustrated  in  figure  3.2.  Note  that  the  critical  point  is  +1  since  all 
feedback  is  considered  positive.  Banklin  [11]  demonstrates  that  the  minimum  distance 
fnnn  any  point  on  the  polar  plot  to  the  critical  point  (the  M  circle  radius),  is  the  reciprocal 
of  the  magnitude  (for  SISO)  or  infinity  norm  (for  MIMO)  of  the  sensitivity.  Since  all 
feedback  is  considered  positive,  the  sensitivity  is  defined  as 

S  =  [I-PKr*  (3.1) 

and  thus  the  M  circle  radius  is 

Mn«i=j^  (3.2) 

where  |S||^,  the  infinity  norm  of  the  sensitivity,  is  the  supremum  over  all  frequency  of  the 
maximum  singular  value  of  L(s). 
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Real(L) 

Figure  3.2.  Polar  Plot  widi  M  Circle  Radius 


PM  =  ±2  arcsin 


(phase  margin  (deg)) 


(3.5) 


The  advantage  of  vector  margins  over  classical  margins  is  that  they  take  into  account 
simultaneous  perturbations  in  both  gain  and  phase,  whereas  the  latter  do  not.  Thus 
classical  margins  may  give  an  overly  t^timistic  prediction  of  the  stability  of  the  system  in 
the  face  of  plant  uncertainty.  Therefore,  vector  margins  are  always  more  ctmservative. 


3-6 


Because  |S|^  can  never  be  less  than  one,  the  vectcv  margins  described  in  equations  (3.3), 

(3.4),  and  (3.5)  have  a  limited  range  of  values  which  are: 

0.0  <  GMy 
—6.0^  GM}  <0.0 

0.0  <  PM  ^  60.0 

This  thesis  will  use  vector  margins  as  a  measure  of  performance  and  a  constraint  for 
the  closed  loop  systems  to  be  developed  QDES.  QDES  allows  easy  constraint  of  diese 
margins  1^  limiting  the  maximum  value  a  specified  SISO  transfer  function  may  take  on.  If 
this  transfer  function  is  the  input  or  output  sensitivity,  then  the  M  circle  radius  and  thus 
the  vector  margins  can  be  ctnistrained  as  desired.  Tire  actual  mechanics  of  diis  process  are 
discussed  in  section  Chapter  VII. 

Summary 

This  chapter  discussed  the  medsures  of  performance  to  be  evaluated  while  using 
QDES  to  explore  the  limits  of  achievable  performance  for  aircraft.  These  measures  were 
divided  into  the  categories  of  aircraft  performance,  time  response  characteristics,  and 
stability  margins.  A  brief  explanation  of  vectcn*  margins  was  also  included. 
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IV.  Classical  Design 


4.1  Introduction 

As  stated  in  Section  (1.5),  a  control  design  for  N2  to  for  a  fighter  aircraft 

was  developed  using  classical  control  methods.  The  purpose  was  two  fold.  First,  the 
results  would  show  what  level  of  perfonnance  in  terms  of  time  response  characteristics 
and  stability  margins  could  be  achieved  manually.  Secondly,  the  classical  design  would 
provide  a  consistency  check  for  the  QDES  results  and  verify  that  the  program  was  being 
used  properly. 

4.2  Axes  System  and  Sign  Convention 

Before  discussing  the  aircraft  model  to  be  used,  it  is  necessary  to  ^wcify  the  axis 
system  and  sign  convention  used  in  the  tmdel.  The  linearired  equations  of  noofitm  for  the 
longitudinal  axis  are  usually  done  in  the  stability  axis  system.  In  this  system,  x^  is  aligned 
with  the  aircraft  velocity  vector  (assuming  no  sideslip),  is  down,  orthogonal  to  X5  and  in 
the  aircraft  plane  of  symmetry,  and  y^  is  directed  out  the  right  wing,  orthogonal  to  x^  and 
Zg.  The  equilibrium  flight  path  angle  (7o)  is  measured  between  the  local  horizon  (Xy  in  the 
vehicle  carried  axis)  and  x^.  The  trim  angle  of  attack  (Otnm)  is  measured  from  x^  to  some 
specified  fixed  reference  on  the  body,  usually  either  the  wing  cord  line  or  the  zero  lift  line. 
These  angles  and  axes  are  illustrated  in  figure  (4.1). 
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vekxity  (y,) 

hx^  hcxizon  X 


Figure  4.1.  Longitudinal  Axis  System 

The  aircraft  pitch  angle  (6),  angle  of  attack  (a),  flight  path  angle  (y),  and  |»tch  rate 
(q)  are  defined  as  positive  in  the  nose  up  direction.  The  elevator  deflection  (5^)  is  positive 
trailing  edge  down  as  shown  in  figure  (4.2). 


Figure  4.2.  Elevator  Sign  Convention 

A  positive  elevator  deflection  then  results  in  a  positive  normal  acceleration  (Nx) 
which  is  nose  down.  This  is  the  normally  accepted  conventicHi.  However,  for  this  thesis, 
positive  N2  will  be  taken  as  nose  up  and  the  aircraft  C  matrix  will  be  modified  to  reflect 
this.  From  a  performance  perspective  this  makes  more  sense  as  most  designers  and  flyers 
are  interested  in  how  many  positive  g's  as  seen  in  the  cockpit  (nose  up)  the  aircraft  can 
pull.  The  ability  of  an  aircraft  to  pushover  mcMPC  than  two  or  three  g’s  is  not  of  great 
operational  importance. 

4.3  Aircraft  Model 

The  aircraft  longitudinal  model  used  m  this  and  the  QDES  (tesigns  is  derived  fiiom 
the  linearized  longitudinal  equations  of  motitHt  at  a  specified  flight  ccmdition  with  it's 
attendant  equilibrium  values.  This  linearization  process  is  thoroughly  developed  in  both 
Etkin  [9]  and  McRuer  [13].  The  assumptitxis  and  simplifications  associated  with  it  were 
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listed  in  Section  (1.3).  The  data  chosen  for  this  design  is  the  A-4D  Skyhawk  at  flight 
condition  5  as  listed  in  McRuer  [13].  Flight  condition  5  parameters  are  in  table  (4.1). 


Altitude  (h) 

(ft) 

Mach  (M) 

Equilibrium 
Velocity  (Vq) 
(ft/sec) 

TrimAOA 

(®trim) 

(deg) 

Eq  Fit  Path 
Angle  (7e) 
(deg) 

15,000 

0.6 

634 

3.4 

0.0 

Table  4.1.  A-^ 

iD  Flight  Condition  5  Parameters 

The  state  space  representation  for  the  longitudinal  axis  is  shown  in  Appendix  (A),  section 
(A.1). 


4.4  Inner  Loop  Pitch  SAS  Theory 

Control  of  an  aircraft  longitudinal  axis  which  exhibits  the  classical  Short  Period  and 

Phugoid  modes  (as  does  the  A-4D)  is  complicated  by  the  very  low  I^ugoid  damping  ratio 
(Cp).  Furthermore,  the  uncompensated  system  does  not  have  a  pole  at  the  origin  of  the  s- 

plane  and  thus  may  exhibit  an  unacceptable  steady  state  error.  Hgure  (4.3)  shows  the  A- 
4D  response  using  a  static  gain  feedback  of  to  a  unit  step  command,  N2(;ind  - 
the  poor  damping  and  aforementioned  steady  state  error.  The  static  gain  feedback  cannot 

correct  the  damping  problem  because,  although  it  drives  the  Phugoid  poles  to  their  zeros, 
decreasing  their  influence  and  increasing  ^p,  it  does  so  at  the  expense  of  Short  Period 

damping  (^sp).  Also,  because  it  adds  no  pole  at  the  origin,  e^  is  quite  poor  (100%).  A 

proportional  plus  integral  (PI)  controller  corrects  the  e^s  problem,  but  suffers  from  the 
same  limitations  in  terms  of  damping.  Thus,  a  method  known  as  Inner  Loop  Pitch 

Stability  Augmentation  System  (SAS)  was  developed.  This  method  uses  an  inner  loop 
static  gain  to  feed  back  pitch  rate  (q)  around  the  plant  to  increase  ^^p.  A  PI  controller  is 

then  placed  in  the  outer  loop  ahead  of  the  plant  to  assure  zero  e^s  and  to  draw  the  locus 
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4  J  Inner  Loop  Pitch  SAS  for  the  A-4D 


(4.5). 


To  choose  the  inner  loop  gain,  stan  with  the  root  locus  of  — —  as  shown  in  figure 

^ecmd 


The  open  loop  poles  of  the  longitudinal  axis  are:  -1.12012  ±  3.50943] 

-0.005512  ±  0.06324] 


-10.0 

Examination  of  the  root  locus  shows  that  the  best  that  can  be  achieved  is  about  0.65. 
This  occurs  at  an  inner  loop  gain  (Kq)  of  -.1883.  The  loop  is  closed  resulting  in  the 
following  inner  loop  closed  loop  poles:  -4.0001  ±  4.74989] 

-.00528  ±  0.05758] 

^.24067 

which  now  become  the  open  loop  poles  of  the  outer  loop. 
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Next,  the  outer  loop  root  locus  of  -  ^  is  examined  to  determine  the  best  location  for  the 

^ecmd 

zero  inserted  by  the  PI  controller  via  the  parameter  Kin2-  Placing  the  zero  at  -8  results  m 

the  root  locus  shown  in  figure  (4.6).  This  placement  pulls  the  locus  further  into  the  LHP. 
Selection  of  the  outer  loop  gain  (K112)  of  0.7  results  in  the  best  combinadtxi  of  Csp  and 

Cp- 

The  closed  loop  poles  are: 

-3.84535  ±  3.42199j 
-2.16893  ±  2.45225j 
-0.007114 
-0.0000003 
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to  a  unit  stq)  input  is  examined. 


The  loop  is  then  closed  and  the  system  tesponse  — — 

^tcmd 

This  time  response  with  its  measures  of  performance  are  ^wn  in  figure  (4.7).  The 
elevator  deflection  (5^)  is  shown  in  figure  (4.8)  while  the  change  in  forward  velocity  (v) 
and  angle  of  attack  (a)  are  shown  in  figure  (4.9). 


Figure  4.7.  A-4D  Inner  Loop  Pitch  SAS  Unit  Step  Re^nse 
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Figure  4.8.  A-4D  Elevator  Deflection  For  1  g  Step  Response 


Figure  4.9.  A-4D  FOTward  Velocity  and  AOA  Change  For  1  g  Step  Response 
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It  is  important  to  remember  that  these  are  perturbation  quantities  away  fimn  the 
equilibrium  values  as  shown  in  table  (4.1).  Thus,  since  the  aircraft  is  in  steady  level  1  g 
flight,  the  actual  N2  the  pilot  feels  is  2  g's.  Likewise,  for  the  angle  of  attack  it  is  necessary 
to  add  the  ottrim  value  to  that  shown  in  figure  (4.9)  to  determine  what  a  the  aircraft 
actually  experiences.  The  same  holds  for  forward  vehxnty  v  change  and  its  equilibrium 
value  (Vq). 

Hgure  (4.10)  shows  the  polar  plot  of  L  for  the  A-4D  with  the  inner  loop  pitch  SAS 
controller.  The  M  circle  radius  of  0.62  results  in  the  following  vector  margins  for  this 
design:  GMy  =  8.4dB 

GMi  =  -4.2  dB 
PM  =  ±36  deg 
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4.6  Summary 


This  chapter  introduced  the  aircraft  model  for  the  longitudinal  axis  and  defined  the 
axis  system  and  sign  ctmvention  used.  The  inner  loop  pitch  SAS  thecny  was  briefly 
explained  and  then  applied  to  the  A-4D  aircraft  N2  to  problrai.  The  results  shown 
provide  a  useful  comparison  with  the  results  obtained  with  QDES  for  the  same  probleno. 
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V.  QDES  Design  for  The  A-4D 


5.1  Introduction 

This  chapter  explores  using  QDES  to  strive  the  Nx  to  control  problem  for 

the  A-4D  as  was  done  classically  in  Chapter  IV.  Widi  the  classical  d^gn,  the  im)cess 
was  performed  and  the  controller  develq)ed  resulting  in  certain  measures  of  performance 
which  were  evaluated  after  the  fact.  The  alrility  to  go  bade  to  say,  achieve  better  stability 
margins  or  a  faster  rise  time,  is  theoretically  possible,  but  would  be  more  of  a  trial  and 
error  process.  With  QDES  the  process  is  reversed.  Starting  with  the  same  (riant  model, 
one  can  specify  several  characteristics  of  the  desired  closed  loop  system  and  then  let 
QDES  attempt  to  develop  the  controller  that  will  deliver  diem.  This  chapter  will  walk 
through  that  process,  ultimately  using  QDES  to  evaluate  the  limits  of  pcrfonnance  for  the 
A-4D  as  defined  in  Chapter  III.  Spedfic  details  of  using  QDES  are  reserved  for  the  step- 
by  step  procedure  outlined  in  Chapter  VII. 

5.2  Model  Development 

The  first  step  in  using  QDES  is  develqiing  the  model.  Through  this  process,  the 
various  input  and  output  vectors  to  the  system  as  described  in  Section  (2.2)  are  identified. 
This  makes  the  necessary  transfer  functions  available  for  QDES  to  use  as  it  seeks  to  find  a 
control  law  to  satisfy  the  s{)ecified  constraints.  For  the  A-4D  example,  the  vectors  are; 

•  Exogenous  Inputs  (w): 

wl:  Nonnal  Acceleration  Command  (Njcmd) 

w2:  Loop  Input  for  StabiliQr  Margin  Constraint/Measurement  (Uin) 

•  Regulated  CXitputs  (z): 

zl:  Angle  of  Attack  (a) 
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z2:  Nonnal  Acceleration  of  the  Aircraft  c.g.  (N2) 

z3:  Pitch  Rale  (q) 

z4:  ElevatOT  Deflection  (5^) 

z3:  Elevator  Deflectiim  Gxmnand  (S^cmd)  stid  Loc^  Output  fcM*  Margins  (Uout) 

•  Sensed  Outputs  (y): 
yl:  Pitch  Rate  (q) 

y2:  Ntmnal  Acceleration  of  the  Aircraft  c.g.  (N2) 
y3:  Normal  Acceleratkm  Conmiand  (N2cjq(1) 

•  Actuator  Inputs  (u): 

ul:  Elevator  Deflection  Command  (S^cmd) 

The  rational  for  these  choices  is  as  follows.  The  transfer  function  between  Uoyt  and  Ujn 
can  be  shown  to  be  the  sensitivity  (S)  as  shown  in  equation  (3.1).  Thus,  making  this 
available  to  QDES  allows  constraint  and  obsovation  of  S,  and  therefore  the  stability 
margins  as  described  in  Section  (3.4).  Outputting  5^  and  a  allows  observation  and 
constraint  of  these  values  which  is  important  for  determining  the  limits  of  performance; 
i.e.,  when  the  stall  a  or  maximum  value  of  5^  is  reached,  the  maximum  attainable  N2  is 
determined.  N2  and  q  were  chosen  for  feedback  as  rate  and  acceleration  senses  tend  to 
provide  better  infcnmation  than  pure  angular  measurements.  N2cix)d  is  also  made  available 
to  the  controller  based  on  DeLaney's  [6]  experiences  with  QDES,  and  also  because  it 
sinply  makes  sense  to  provide  this  known  signal  to  the  controller.  The  choice  of  N,nniH 
as  an  input  and  N2  as  an  output  is  obvious. 

Figure  (5.1)  shows  the  block  diagram  representation  of  the  aircraft  plant  and 
controller  with  these  input  and  output  vectors.  Note  that  this  diagram  represents  the  same 
relaticmship  between  the  input  and  output  vectors  and  the  plant  components  as  the 
standard  feedback  diagram  of  figure  (2.1).  Howevo-,  diis  representation  makes 
construction  of  the  nine  matrices  used  in  equations  (2.2),  (2.3),  and  (2.4)  much  more 
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straightf(»ward.  In  figure  (S.l),  G  is  tlM  A-4D  state  q>8ce  widiout  die  actuator  as  shown 
in  Appendix  A,  Section  (A.1.1),  and  is  die  actuator  state  ^Mce,  shown  in  Section 
(A.1.2).  K|  and  K2  together  represent  the  single  K(s)  in  figure  (2.1). 


In  equations  (S.2)  and  (S.3),  Cm  and  refer  to  the  "nz"  row  the  C  and  D  matiices  of 
G  in  figure  (S.l). 

5.3  Plant  Discretization 

As  described  in  Section  (2.7),  the  plant  to  be  controlled  QDES  must  be 
discretized.  For  this  thesis,  diis  was  easily  accomplished  by  use  to  the  "c2dm"  ctmunand 
in  MATLAB.  The  only  decision  that  needs  to  be  made  is  the  discretization  rate.  This  is 
driven  by  the  highest  frequency  mode  in  the  plant  In  Qiapter  IV,  the  A-4D  plant  poles 
were  identified.  The  conesponding  firequencies  are:  .06  rad/sec  (.0095  Hz) 

3.7  rad/sec  (.5888  Hz) 

10  rad/sec  (1.592  Hz) 

The  highest  frequency  ccxnes  from  the  actuate  pole.  Houpis  [12]  recommends  sampling 
at  least  twice  as  fast  as  the  highest  Oequency  (Shannon's  Sampling  Theorem).  However,  it 
has  been  suggested  that  sampling  at  flve  times  the  highest  frequency  is  actually  better  [14]. 
Thus,  for  the  A-4D,  a  sampling  rate  of  8  Hz  was  used.  Certainly,  a  higher  sampling  rate 
can  always  be  used,  but  there  is  litde  utility  in  doing  so  as  is  discussed  in  Chapter  VII. 
The  plant  that  was  discretized  is  the  A-4D  as  listed  in  Appendix  (A),  Section  (A.  1.3), 
which  includes  the  actuator  dynamics,  and  has  been  scaled  to  the  proper  units  with  N2 
defined  as  positive  up. 

5.4  Generating  the  Impulse  Response  Matrices 

As  described  in  Section  (2.7),  QDES  starts  with  the  in:q>ulse  responses  of  the  three 
matrices  (Tl,  T2,  and  T3)  which  describe  the  closed  loop  plant  with  a  nominal  stabilizing 
controller  as  represented  in  equation  (2.29).  To  generate  these  for  the  A-4D  example,  a 
MATLAB  script  file  "nzprp.m"  was  written  as  shown  in  Appendix  (C),  Section  (C.1). 
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The  file  is  heavily  commenied,  identifying  each  step  of  Um  process  with  instructions  as  to 
what  infcvniation  is  to  be  inserted  for  a  generic  aircraft  plant  model.  The  following  is  a 
walk  through  this  file  with  explanations  and  rationale  for  the  paianKtm  and  processes 
used. 

The  original  plant  model  as  listed  in  Section  (A.1.1)  is  entered.  The  actuator 
dynamics  are  augmented  to  the  plant  and  N2  is  changed  to  the  positive  up  convention 
resulting  in  the  matrices  in  Sectitm  (A.  1.2).  Fiiuilly,  the  units  are  converted  resulting  in  the 
state  space  shown  in  Section  (A.  1.3). 

The  scale  factor  is  used  to  change  the  value  of  N2cm(i.  Since  QDES  ctmstrains  the 
tioK  domain  characteristics  in  terms  of  the  unit  step  response,  the  sysrem  needs  to  be  told 
that  a  unit  step  input  is  really  another  value,  say  4  g's.  By  dividing  the  system  C  matrix  by 
the  scale  factor,  QDES  designs  a  closed  loq}  system  as  if  the  commanded  value  was  the 
scale  factor  times  the  unit  step.  The  various  constraints  such  as  a  and  5^  are  properly 
interpreted  while  the  internal  states  reflect  their  actual  values  for  the  non-unity  input 

The  plant  is  discretized  using  the  "c2dm"  command  as  previmisly  discussed.  Note 
that  the  command  asks  for  the  sample  time  (Ts)  which  is  the  inverse  of  the  sample  rate. 
Also,  the  type  of  discretization  must  be  specified.  In  this  case,  the  Zero  Order  Hold  ('zdi') 
was  used. 

The  nine  matrices  as  developed  in  Section  (5.2)  are  entered.  Notice  that  wherever 
an  entry  is  not  either  a  one  or  a  zero,  the  data  fiom  the  discretized  plant  is  inserted.  Also, 
die  matrices  must  be  labeled  as  listed  so  as  to  be  properly  read  the  MATLAB  script  file 
"obc2ttt.m". 

As  previously  discussed,  a  nominal  stabilizing  controller  must  be  designed  to  form 
the  nominal  closed  loop  plant.  There  are  several  ways  to  do  this.  The  method  en^loyed 
here  involves  choosing  a  full  state  feedback  gain  matrix  (K)  using  the  MATLAB  "dlqr" 
command  and  the  Kalman  filter  gain  matrix  (L)  using  the  command  "dlqe".  These 
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matrices  are  then  comlnned  with  the  plant  infonnatkm  entered  in  the  nine  matrices  using  a 
parameterization  similar  to  that  listed  in  equations  (2.33),  (2.34),  and  (2.3S)  to  ftnm  the 
nominal  closed  lot^  system.  To  use  "dlqr",  two  weighting  matrices  must  be  chosen.  The 
first  is  the  state  weighting  matrix  Ql.  The  method  employed  is  the  same  as  used  by 
DeLaney  [6]  whereby  the  states  are  weighted  by  the  inverse  of  the  square  of  the  maximum 
value  the  state  is  allowed  to  take  on.  Fot  exanq)le,  the  A'4D's  elevator  is  considered  to  be 
limited  to  ±  25  deg  .  Therefore,  the  5^  state  is  weighted  by  0.(X)16.  The  state  weighting 
matrix  for  this  example  is: 


0.00044 

0 

0 

0 

0  ■ 

0 

0.0016 

0 

0 

0 

0 

0 

0.00111 

0 

0 

0 

0 

0 

0.00111 

0 

0 

0 

0 

0 

0.0016  _ 

Note  that  this  matrix  must  be  positive  semidefinite  meaning  that  its  eigenvalues  are  ^  0. 
The  other  matrix  to  be  chosen  is  the  control  weight.  This  value  is  chosen  somewhat 
arbitrarily,  and  in  the  absence  of  any  concrete  guidance,  DeLaney's  number  of  0.1  was 
used.  Thus,  the  control  weighting  matrix  is; 

Rl  =  [0.1] 

This  matrix  must  be  positive  definite  (eigenvalues  >  0)  to  prevent  the  use  of  infinite 
control  power.  To  use  "dlqe"  more  weightings  must  be  chosen.  The  first  is  G  which 
determines  what  states  could  be  affected  by  process  noise.  For  an  aircraft,  this  might  be  a 
wind  gust.  A  value  of  unity  is  put  on  every  state  that  could  be  affected  by  such  a  gust,  and 
a  zero  on  each  that  could  not.  G  for  this  design  is: 
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1 


G  = 


1 

1 

1 

0 


The  other  two  weights  to  be  specified  are  Q2  and  R2.  These  represent  the  expected  value 
of  the  state  and  measurement  noises  respectively.  Again,  using  DeLaney's  numbo^: 

Q2  =  [0.1] 


R2  = 


0.001 

0 

0 


0 

0.001 

0 


0 

0 

0.001 


The  choice  of  weightings  is  somewhat  arbitrary.  In  theory,  as  long  as  the  mminal  system 
is  stable,  its  actual  performance  characteristics  are  not  that  important.  QDES  should 
arrive  at  the  same  closed  loop  system  (if  it  can  find  a  feasible  solution  at  all),  although  it 
may  take  more  or  less  time  depending  on  how  close  the  nominal  system  is  to  the  desired 
system.  As  will  be  seen,  this  was  not  necessarily  the  case. 

Once  the  weightings  are  chosen  to  derive  K  and  L,  two  other  pieces  of  information 
need  to  be  made  available  to  "obc2ttt.m"  before  it  can  be  called.  The  first  is  the  source  file 
name.  This  name  must  match  the  QDES  source  file  name  and  identifies  the  in^ulse 
response  matrices  Ctl.  _t2,  _t3)  plus  the  "_tttmat"  file  with  the  current  problem.  For  this 
example,  A4D7NZ  was  used.  Also,  the  number  of  samples  must  be  specified.  This  is  the 
plant  discretization  rate  (Hz)  times  the  time  interval  (sec)  over  which  you  desire  QDES  to 
solve  the  problem.  The  time  period  chosen  for  the  design  was  five  seconds.  The  rationale 
for  this  choice  is  that  the  purpose  is  to  examine  the  maximum  g  attainable  in  addition  to 
whether  a  controller  can  be  designed  to  work  at  this  level.  Thus,  compliance  with  the  time 
and  margin  ctmstraints  is  as  important  as  finding  when  the  5^  and  a  limits  are  reached. 
Additionally,  since  control  of  thrust  is  not  included  in  the  model,  examination  of  sustained 
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g  capability  is  not  reaUy  possible.  Lastly,  sustaining  even  moderate  g  for  an  extended 
period  of  time  would  result  in  a  change  m  pitch  angle  (6)  outside  that  valid  for  the  small 
angle  approximation  used  in  the  linearization  process  alluded  to  in  Section  (1.3).  Thus, 
for  this  example,  a  rate  of  8  Hz  for  S  sec  yields  40  samples. 

Finally,  the  MATLAB  script  file  ”obc2ttt.m"  is  called.  This  file  actually  calls  the 
function  "obc2ttt_fn.m".  Note  that  to  use  ''obc2ttt.m",  the  plant  must  be  strictly  prq>er, 
that  is  pydu  must  equal  zero.  Given  the  state  space  representations  used  in  this  thesis  with 
the  actuator  dynamics  augmented,  this  will  always  be  the  case.  These  files  are  siq)plied 
with  the  QDES  tape  as  described  in  Erickson  [8]  as  are  the  others  necessary  to  carry  out 
the  process  described  in  this  thesis.  The  qrecifics  of  what  files  are  needed  and  how  to  use 
them  is  described  in  Chapter  VII.  When  the  file  ''a4d7nzprp.m"  was  mn,  the  A4D7NZ_tl, 
A4D7NZ_t2,  A4D7NZ_t3,  and  A4D7NZ_tttmat  files  where  output  for  later  use  by 
QDES. 

5.5  Writing  the  QDES  Source  File 

Once  the  inqrulse  response  matrices  have  been  generated,  the  QDES  source  file 
must  be  written  to  specify  the  problem  so  QDES  may  attempt  a  solution.  The  source  file 
is  written  in  "C"  and  close  attention  to  the  examples  in  the  QDES  User  Craimands  [4]  is  a 
must  for  avoiding  syntax  errors.  The  source  file  is  divided  into  four  parts  which  will  be 
discussed  separately.  The  source  file  for  this  A-4D  exaiiq>le  is  shown  in  Appendix  (C), 
Section  (C.2). 

5.5.1  The  Preamble  In  the  preamble,  the  exogenous  inputs,  regulated  outputs, 
and  sensed  outputs  as  listed  in  Section  (5.2)  are  defined  by  a  symbol  of  the  user's  choice. 
They  are  also  numbered,  and  this  numbering  must  reflect  the  order  of  the  variables  listed 
in  Section  (5.2)  and  as  used  in  equations  (5.1),  (5.2),  and  (5.3). 
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5J.2  The  Declarations  Section  This  section  tells  QDES  how  many  of  each  of  the 
variables  listed  in  the  preamble  as  well  as  the  number  of  actuators  are  present  in  the 
problem.  Also  input  is  the  number  of  samples  from  the  "nzprp.m"  file  and  the  number  of 
FIR  taps.  The  number  of  taps  determines  the  number  of  terms  in  the  finite  difference 
equations  which  are  the  discrete  equivalents  of  the  continuous  state-variable  diffeieiitial 
equation  as  shown  in  Chapter  II.  More  taps  equals  more  information  that  QDES  has 
available  to  solve  the  control  problem.  There  are  restrictions  to  the  number  of  taps  used 
however.  Boyd  [4]  recommends  limiting  the  product  of  the  number  of  actuators,  number 
of  sensed  outputs,  and  the  number  of  taps  to  less  than  200  due  to  limitatitnis  of  the  convex 
problem  solver  used.  For  the  problem  set  up  as  used  in  this  thesis,  that  number  is  66. 
Furthermore,  £)eLaney  [6]  states  that  the  number  of  taps  should  be  strictly  less  than  the 
number  of  samples,  or  unconstrained  dynamics  may  appear  in  the  Q  portion  of  the 
compensator.  In  a  few  attempts  where  this  restriction  was  violated,  the  author  did  not 
experience  this  phenomenon.  Lastly,  the  size  of  the  QDES  controller  is  equal  to  the  size 
of  the  nominal  controller  plus  the  number  of  taps.  Thus,  if  the  size  of  the  controller 
required  is  one  of  the  elements  of  the  design  investigation,  this  value  could  be  adjusted 
accordingly.  More  discussion  of  the  choice  of  number  of  taps  follows  in  Chapter  VII.  In 
this  example,  35  taps  were  used.  The  last  entry  in  this  section  are  the  names  of  the 
impulse  response  matrices  that  were  created  with  "nzprp.m"  using  "obc2ttt.m". 

5.5.3  The  Objective  Section  This  section  tells  QDES  to  minimize  one  of  several 
available  relationships  between  the  inputs  and  outputs,  or  a  time  response  characteristic  as 
described  in  [4].  Use  of  objective  functions  was  tried,  but  not  used  for  the  final  results  of 
this  thesis  for  several  reasons.  First,  based  on  his  experience  with  QDES,  DeLaney 
recommended  not  using  them  due  to  the  time  it  takes  for  the  program  to  run.  The  author 
found  this  to  be  a  problem  as  well.  Secondly,  as  the  purpose  here  is  to  evaluate  the 
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maximum  N2  attainable  by  the  aircraft  given  a  set  of  omstraints,  there  was  no  real  need  to 
attempt  to  minimize  any  of  these  constraints,  and  thus  using  an  objective  functicm  only 
served  to  slow  the  process  down. 

5J.4  The  Constraints  Section 

The  constraint  section  allows  the  user  to  impose  limits  on  the  closed  loop  system 
to  be  designed.  These  take  the  form  of  time  response  characteristics  and  various 
input/output  relationships,  again  as  described  in  [4].  For  this  thesis,  three  different  time 
response  constraints  were  used  along  with  a  constraint  on  the  magnitude  of  the  transfer 
function  representing  the  sensitivity.  The  complete  list  of  the  constraints  used  is  found  in 
the  source  file  in  Section  (C.2).  A  thorough  discussion  of  the  mechanics  of  writing  the 
constraint  section  is  found  in  Chapter  VII.  For  purposes  of  this  chapter,  the  time  respmse 
constraints  for  this  example  were:  T,  <  1.5  sec 

Mp<10% 

Us  <  5  % 

Tji  <  2.5  sec  (3  %) 

®cin(i  <  3  % 

These  constraints  form  an  envelope  within  which  the  time  response  is  to  remain.  This  is 
seen  on  the  time  response  graphs  that  follow.  These  were  chosen  based  on  the  results 
from  the  classical  design,  but  relaxed  slightly  especiaUy  for  e^^j.  Experience  with  QDES 
showed  that  making  e^^d  needlessly  restrictive  often  resulted  in  no  solution.  The  elevator 
deflection  was  constrained  to  ±  25  deg  and  angle  of  attack  as  follows:  -5  deg  <  a  <  18 
deg.  Note  these  are  not  necessary  the  actual  A-4D  limits  but  are  representative  of  aircraft 
of  this  type.  Also,  it  must  be  remembered  that  the  values  derived  from  the  linearized 
model  and  thus  from  evaluation  of  QDES  results  are  perturbation  values.  Thus,  if  the 
aircraft  has  an  angle  of  attack  limit  of  20  deg  aixi  the  flight  condition  atri„,  is  1.5  deg,  the 
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upper  constraint  would  be  given  to  (^ES  as  18.5  deg.  Lastly,  the  vector  margins  [-4,7] 
dB  and  [±3S]  deg  were  imposed  by  constraining  |S|„  <  1.6628.  Actually,  if  the  sensitivity 

constraint  is  met,  the  vector  gain  margins  will  be  better  than  those  shown.  This  possiUy 
confusing  relationship  is  fuUy  explained  in  Cht4}ter  VII.  The  choice  of  margin  constraints 
was  again  driven  by  the  results  of  the  classical  design  and  the  fact  that  the  lowest 
acceptable  phase  margin  is  generally  considered  to  be  30  deg  [11]. 

5.6  Running  QDES 

Once  the  impulse  response  matrices  are  developed  and  the  QDES  source  file  is 
written,  the  actual  running  of  QDES  is  a  sinq)le  matter.  The  q>ecific  mechanics  of  doing 
so  is  reserved  for  Chapter  VII.  When  the  QDES  mn  is  con^lete,  QDES  fmms  the 
*'sourcefile.q"  and  "sourcefile.out"  files  which  are  subsequently  used  to  form  the  QDES 
controller.  If  QDES  cannot  find  a  controller  that  meets  the  constraints,  it  will  generate  a 
message  to  that  affect  and  output  ”sourcefile.lastq’'  and  "sourcefile.out”  files.  The  ’’.out” 
file  can  be  examined  to  see  what  constraints  QDES  could  not  meet.  Run  times  for  the 
examples  in  this  chapter  were  on  the  order  of  one  minute  or  less,  which  works  weU  with 
the  iterative  process  necessary  to  determine  the  N2  limits  of  performance.  Other 
combinations  of  numbers  of  samples,  number  of  taps,  and  sample  rates  which  were  run  to 
explore  QDES'  capability  resulted  in  significantly  Icmger  run  times. 

5.7  Forming  the  QDES  Controller  and  Closed  Loop  System 

The  QDES  controller  and  closed  loop  system  can  be  easily  formed  via  the 
MATLAB  script  files  "qobc2comp.m'’  and  "gethzw.m"  supplied  with  the  QDES  tape  and 
described  in  [8].  To  employ  these,  a  MATLAB  script  file  ’’nzpstm"  was  written  which  is 
found  in  Appendix  (C),  Section  (C.3).  Again  this  file  is  heavily  commented  and  explained 
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in  detail  in  Chapter  VII.  It  does  require  entering  some  of  the  same  parameters  as  those 
entered  in  the  "nzprp.m”  file.  The  file  also  simulates  the  respmise  of  the  (^ES  closed 
loop  system  (H2w)  to  a  unit  step  and  outputs  graphs  of  the  N2  time  response,  the  5^  and  a 
time  histories,  the  Bode  magnitude  plot  of  the  sensitivity,  and  the  Bode  magnitude  of  N2 
to  Njcmd  closed  loop  system. 

5.8  A-4D  QDES  Design  Results 

This  secticm  examines  the  results  of  using  QDES  to  determine  the  N2  limits  of 
performance  subject  to  the  constraints  listed  in  Section  (5.4.4). 

5.8.1  A-4D  QDES  Design,  1  g  Response  Figure  (5.2)  shows  the  time  response  of 
the  A-4D  N2  to  N2cind  for  1  g.  Notice  that  this  is  not  the  picture  of  an  ideal  step  response, 
but  that  it  does  remain  within  the  time  response  constraint  envelope.  QDES  only 
guarantees  that  the  closed  loop  system  meets  the  constraints,  not  that  the  controller 
developed  produces  a  nice  response.  This  is  also  consistent  with  the  adnKMiitions  of  the 
program's  author,  Dr.  Boyd,  to  the  affect  that  QDES  was  never  intended  to  be  used  to 
design  controllers.  Also,  notice  at  the  five  second  point,  the  response  is  just  outside  the 
constraint  envelope.  This  is  due  to  the  wording  of  the  constraint  which  requires  the 
response  to  meet  the  constraints  up  to  the  n_sample-l  point. 
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Figure  (5.3)  shows  the  time  histories  of  the  elevator  deflection  and  the  angle  of  attack. 
From  the  graphs,  one  can  estimate  that  the  g  response  wiU  probatdy  be  a-limited  at  about 
4.5  g's  when  the  a  response  hits  its  limits  of  18  deg.  QDES  will  verify  whether  a 
controller  can  be  designed  at  this  level  and  if  any  changes  in  the  response  can  be  made  to 
achieve  a  higher  level  of  performance. 


5-13 


The  top  half  of  figure  (5.4)  shows  the  Bode  magnitude  of  the  input  sensitivity  in  absdute 
magnitude,  not  in  dB.  The  maximum  value  attained  on  this  plot  is  |S|„  fiom  which  the 

following  vector  margins  for  this  design  were  calculated: 

GM„=-2.6dB 
GMi=+3.8dB 
PM  =  ±20  deg 

Notice  that  the  maximum  value  of  the  input  sensitivity  is  greater  than  the  value  (1.6628) 
specified  in  the  constraint  section  of  the  QDES  source  file  even  though  QDES  indicated 
that  it  had  met  all  the  constraints  by  returning  a  solution.  Chapter  VIII  addresses  diis 
l^enomenon  and  a  theory  as  to  why  it  occurs.  The  bottom  half  of  figure  (5.4)  shows  the 
Bode  magnitude  (dB)  of  the  closed  loop  N2  to  N2(.|nd  transfer  fiinction.  This  is  included 
to  examine  the  closed  loop  bandwidth  (CLBW).  Whoi  dealing  with  a  discrete  system, 
Anderstm  [1]  suggests  that  one  needs  to  look  at  the  CLBW  of  the  final  design  to  ensure 
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that  the  ratio  of  the  sanq^g  rate  to  the  CXBW  is  between  about  4  and  20.  For  a 
sampling  rate  of  8  Hz  (SO  rad/sec),  the  CLBW  would  have  to  exceed  q^troximately  12 
rad^  before  the  lower  end  of  diis  range  was  neared.  As  seen  in  the  plot,  the  (XBW  is 
approximately  2.S  rad/sec  (-3dB  down  point).  Thus,  the  sample  rate  of  8  Hz  is  still  valid. 


5.8.2  A-4D  QDES  Design,  3  g  Response  Hgures  (5.5),  (5.6),  and  (5.7)  show  the 

same  results  for  the  3  g  response  of  the  QDES  design.  The  cmly  marked  difference 

between  these  results  and  those  for  1  g  occurs  with  the  input  sensitivity.  In  diis  case  it  is 

omstrained  to  be  less  than  1.6628  resulting  in  the  following  vector  margins; 

GM„=-4.3dB 
GMi=+8.7dB 
PM  =±37  deg 

Again,  Chapter  VUI  will  examine  this  occurrence. 
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5.8.3  A~4D  QDES  Design,  45  g  Response  Figures  (5.8),  (5.9),  and  (5.10)  show 

the  results  of  the  4.5  g  QDES  design.  Notice  an  expansion  of  the  a  plot  shows  it  is  right 

against  the  18  deg  limit  One  might  expect  that  the  maximum  N2  attainable  for  these 

constraints  has  been  reached.  According  to  QDES,  this  is  the  case  as  4.6  g's  was  tried  and 

QDES  could  not  return  a  solution.  The  stability  margins  for  this  design  are: 

GMu=-4.1dB 
GMi=+8.1dB 
PM  =±35  deg 
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alpha  (deg)  Normal  g  (gs) 


5.8.4  A-4D  QDES  Design,  4.6  g  Response  In  an  attempt  to  explore  how  QDES 
would  modify  its  closed  loop  system  to  changes  in  the  constraints,  T,  and  Tg]  were  relaxed 
by  0.5  sec  and  a  4.6  g  problem  was  tried.  Figures  (5.11),  (5.12),  and  (5.13)  show  the 
results.  Note  how  the  system  displays  some  nonminimum  phase  activity  and  has  reduced 
Mp  such  that  a  does  not  spike  above  18  degs  initiaUy  thus  allowing  a  4.6  g  response.  The 
gain  and  phase  margins  are  identical  to  the  4.5  g  example. 
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a^a(deg)  _  deltae(deg)  Normal  g(gs) 
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Figure  5. 13.  A-4D  4.6  g  Input  Sensitivity  and  CLBW 


5.9  A-4D  Short  Period  Approximation 

Since  the  dme  period  over  which  the  response  of  the  aircraft  is  being  examined  is 
relatively  short,  it  is  of  interest  to  see  what  QDES  will  do  using  the  Short  Period 
approximation.  For  the  Short  Period  approximation,  the  I^ugoid  mode  is  removed  since 
its  response  is  manifested  at  relatively  long  time  periods  after  the  q>plication  of  the  input 
During  the  Phugoid  part  of  the  response,  the  forward  velocity  (v)  changes  very  little. 
Therefore,  to  form  the  Short  Period  approximadem,  the  v  row  and  column  are  renKwed 
fiom  the  A  matrix  and  the  v  row  is  removed  from  the  B  matrix  in  the  A-4D  model  in 
Section  (A.  1.1).  When  this  is  done,  the  pitch  rate  (6)  column  of  the  A  matrix  is  all  zeros, 
and  the  state  space  realization  is  non-minimal;  i.e.,  the  6  state  is  not  needed.  Therefore  it 
is  removed  as  well.  The  resulting  state  space  realization  with  the  actuator  dynamics  and 
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scaled  to  the  proper  units  is  listed  in  Appendix  (A),  Section  (A.3).  The  state-variable 
differential  equations  for  the  A-4D  Short  Period  approximation  are: 


All  processes  and  procedures  for  using  QDES  are  the  same  as  presented  for  the  full  A-4D 
model  except  that  the  Short  Period  plant  was  discretized  at  20  Hz  and  50  FIR  taps  were 
used.  This  provided  a  smoother  response  than  that  seen  in  the  full  A-4D  model.  Also, 
QDES  was  constrained  to  provide  zero  command  error.  This  constraint  was  tried  with  the 
full  model,  but  it  did  not  improve  the  response  as  it  did  with  the  Short  Period 
approximation  and  often  it  caused  QDES  not  to  return  a  solution. 

5.9.1  A-4D  Short  Period  QDES  Design,  I  g  Response  The  1  g  response  for  the 
A-4D  Short  Period  QDES  design  is  shown  in  figures  (5.14),  (5.15),  and  (5.16).  The 
differences  from  the  full  model  include  a  smoother  resnonse  that  more  closely  and  steadily 
tracks  commanded  value  and  the  sensitivity  which  lacks  the  low  frequency  spike  shown  in 
the  full  model.  This  is  due  to  the  absence  of  the  low  frequency  Phugoid  pole.  The  vector 
margins  are: 
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1  2  3  4  5 

Time  (sec) 

Figure  5.14.  A-4D  Short  Period  QDES  Design,  1  g  Response 


Magnitude  (dB)  Magnitude  (not  dB) 


Figure  5.15.  A-4D  Short  Period  Elevator  and  AOA  1  g  Responses 


Figure  5.16.  A-4D  ShOTt  Period  1  g  Input  Sensitivity  and  CLBW 
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5.9.2  A-4D  Short  Period  QDES  Design,  3  g  Response  Figures  (5.17),  (5.18), 
and  (5.19)  show  the  A-4D  Short  Period  QDES  design  for  a  3  g  input.  The  vector  margins 


Magnitude  (dB) 

I  I  I 


0  1  2  3  4  5 

Time  (sec) 

Figure  5.18.  A-4D  Short  Period  Elevator  and  AOA  3  g  Responses 


io’  10°  io’ 

Frequency  (rad^sec) 

Figure  5.19.  A-4D  Short  Period  3  g  Input  Sensitivity  and  CLBW 
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5.9.3  A-4D  Short  Period  QDES  Design,  4.6  g  Response  Figures  (S.20),  (5.21), 
and  (5.22)  the  A-4D  Short  Period  QDES  design  for  a  4.6  g  conunamL  Notice  again 
that  the  response  appears  to  be  very  close  to  the  a  limit  of  18  deg.  In  fact,  it  is  the  limit 
for  these  cemstraints,  as  4.7  g's  was  tried  and  QDES  could  not  provide  a  solution.  The 
vector  margins  are: 

GMu=-4.1dB 
GMi=+8.2dB 
PM  =  ±35  deg 


Figure  5.20.  A-4D  Short  Period  QDES  Design,  4.6  g  Response 
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Figure  5.21 .  A-4D  Short  Period  Elevator  and  AOA  4.6  g  Responses 


10’  10®  10’ 

Frequency  (rad/sec) 

Figure  5.22.  A-4D  Short  Period  4.6  g  Input  Sensitivity  and  CLBW 
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5.9.4  A-4D  Short  Period  QDES  Design,  4.7  g  Response  As  was  dc»c  for  the  full 
model  A-4D  case,  the  T,  and  T,]  constraints  were  relaxed  by  O.S  sec  and  a  QDES  solutitm 
was  again  sought.  As  no  solution  was  found,  the  constraints  were  relaxed  another  O.S 
seconds  with  the  same  result.  Hnally,  the  e^^  =  0  constraint  was  removed  and  a 

solution  found  which  is  shown  in  figures  (5.23),  (5.24),  and  (5.25).  Notice  how  the 

response  hugs  the  lower  3  %  error  envelope  and  how  QDES  has  reduced  Mp  to  meet  the 

constraint  The  vectOT  margins  for  this  design  are: 

GMu=^.ldB 
GMi=-»-8.2dB 
PM  =  ±35  deg 
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5.10  Sununary 


This  chapter  demonstrated  the  use  of  QDES  in  evaluating  the  limiK  d  peifonnance 
for  Nj  for  the  A-4D  aircraft.  It  stepped  through  the  process  of  model  deveh^ment  and 
implementation  with  QDES  through  all  the  associated  MATl-AB  script  files  needed  to 
prepare  the  problem  for  QDES  and  examine  its  solutitm.  A  more  detailed  step-by-step 
procedure  is  provided  in  Chapter  VII. 

The  results  were  fairly  straightforward  with  a  few  interesting  results.  In  the  case  of 
the  full  nmdel,  the  time  response,  although  not  as  smooth  as  in  the  classical  design, 
performed  as  well  in  terms  of  the  time  domain  measures  of  performance  except  for  the 
command  error  .  This  appears  to  be  caused  by  the  Phugoid  nude  not  being  sufficiently 
controlled  by  QDES.  This  is  also  apparent  in  the  lack  of  constraint  of  the  sensitivity  and 
thus  the  vector  margins  for  the  1  g  case.  However,  for  the  higher  g  runs,  the  full  model 
had  margins  within  the  specified  constraints.  The  Short  Period  approximation,  having 
gotten  rid  of  the  Hiugoid  mode,  performed  quite  well  for  all  runs  and  even  managed  a 
slightly  higher  g  level  before  hitting  the  a  limiL  Chapter  VIII  will  examine  the  sensitivity 
result  more  closely. 

For  a  relatively  easy-to-control  aircraft  like  the  A-4D,  the  value  of  using  QDES  may 
not  be  quite  overwhelmingly  apparent.  However,  the  ability  to  re-design  the  system  to 
meet  ever  varying  constraints  for  the  evaluation  of  limits  of  performance  lends  quite  a  bit 
of  flexibility  to  the  design  and  evaluation  process.  Countless  "what  ife"  can  be  examined 
in  a  relatively  easy  manner  compared  to  re-designing  and  tweaking  a  classically  designed 
controller. 
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VI .  QDES  Design  for  the  F-I6 


6.1  Introduction 

This  chapter  explwes  using  (^ES  to  solve  the  N2  to  Njcmd  problem  for  the  F-16  as 

was  done  for  the  A-4D  in  Chapter  V.  Except  where  noted,  the  problem  set  up  and 

constraints  are  the  same  as  those  used  for  the  A-4D.  The  F-lb  Imgitudinal  model  at  sea 

level  and  Mach  0.6  has  the  following  poles; 

-20.0 
-4.3542 
+1.9055 
-0.00861. 0708i 

Note  that  the  RHP  poles  makes  this  aircraft  statically  unstable  and  therefore  much  more 
difficult  to  control.  It  was  of  interest  to  see  what  QDES’  capabilities  were  with  this  model 
as  opposed  to  the  A-4D  which  is  statically  stable.  Appendix  (B)  contains  the  state  q)ace 
representation  which  is  used  in  the  problem  set  up  as  depicted  in  equations  (5.1),  (5.2), 
and  (5.3). 

6 .2  Differences  From  the  A  -4D  Problem 

The  highest  frequency  pole  for  the  F-16  is  the  actuator  at  20  rad/sec.  Note  this  is 
twice  as  fast  as  that  for  the  A-4D,  therefore  the  minimum  sanq)le  rate  must  also  be  twice 
as  fast,  or  16  Hz.  When  run  for  5  sec  this  equates  to  twice  the  number  of  sample  times 
which  does  not  merely  double  the  size  of  the  problem.  In  fact,  it  doubles  the  number  of 
each  of  the  various  constraints  over  their  particular  time  span.  In  the  A-4D  exan:q)le,  the 
job  sizes  were  i^proximately  6  megabytes.  With  the  F-16  using  50  taps  the  job  size  was 
almost  85  megabytes,  which  caused  the  computations  to  slow  to  a  virtual  standstill.  In 
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fact,  above  iQ)proxiinately  4S  megabytes,  the  coaq>uier  entered  a  motte  whereby  it  had  to 
wait  for  sufficient  disk  space  to  become  available  before  perftmning  ctnnputations. 
DeLaney  recommended  never  allowing  the  problem  to  exceed  100  megabytes  when  not 
using  any  objective  functicms.  However,  to  get  reasonable  run  times  (30  minutes  or  less) 
which  would  allow  a  designer  to  quickly  evaluate  his  or  her  preliminary  design,  it  was 
found  the  job  size  neetted  to  be  restricted  to  apivoximately  43  megabytes.  To  accomplish 
this,  either  the  number  of  samples  or  the  number  of  taps  needed  to  be  reduced.  Through 
much  trial  and  error,  it  was  found  that  reducing  the  number  of  taps  excessively  precluded 
a  solution  where  one  was  known  to  exist  from  previous  trials.  Thus,  since  the  sample  rate 
needed  to  remain  at  its  miiiimum  value  of  16  Hz,  the  run  time  was  reduced  to  reduce  the 
number  of  samples.  This  resulted  in  the  F-16  problems  being  nm  at  16  Hz  for  3.5  sec  with 
45  FIR  Taps.  This  resulted  in  job  sizes  of  around  42  megabytes  which  for  the  F-16  full 
model  ran  for  5  to  15  minutes.  Since  this  investigation  is  concerned  with  finding  the 
maximum  attainable  within  the  time  response  and  margin  constraints  and  not  maximum 
sustained  g,  this  change  does  not  impact  the  value  of  the  experiment 

The  other  difference  from  the  A-4D  exanqjle  had  to  do  with  the  constraint  of  the 
sensitivity.  With  the  F-16,  QDES  would  not  return  a  solution  if  S  was  constrained  below 
a  value  of  25.  Thus  all  runs  were  made  with  this  very  liberal  constraint  in  place. 

6.3  F-16  QDES  Design  Results 

The  results  from  the  F-16  design  runs  are  presented  in  the  same  manner  as  for  the 
A-4D.  Significant  differences  will  be  noted  as  they  are  presented  with  a  summary  at  the 
end  of  the  chapter. 

6.3.1  F-16  QDES  Design,  1  g  Response  Rgures  (6.1),  (6.2),  and  (6.3)  show  the 

results  of  the  QDES  design  for  a  1  g  command  in  the  same  format  as  that  used  for  the  A- 


6-2 


4D.  Note  the  N2  response  is  similar  to  the  A-4D  as  are  the  5^  and  a  resptHises,  although 

the  F-16  at  this  flight  cmdidon  uses  less  elevator  and  ^ows  less  increase  in  AOA. 

Imposing  the  same  constraints  as  on  the  A-4D,  one  expects  to  be  able  to  attain  a  higher  g 

level  than  the  A-4D.  In  fact,  it  appears  that  the  F-16  will  also  be  a-limited,  but  at 

approximately  9.8  g's.  Examinadon  of  the  sensidvity  plot  shows  the  same  behavior  as  the 

A-4D  at  Ig.  There  is  a  spike  at  the  Phugoid  finequency  of  0.07  rad/sec  which  exceeds  the 

constraint  value  even  though  QDES  claimed  to  provide  a  soludrai  that  met  all  constraints. 

The  rest  of  the  plot  above  the  Hiugoid  spike  is  constrained  below  25  as  requested.  This 

spike  at  the  Phugoid  frequency  results  in  predictaUy  poor  margins  of: 

GMu=-0.24dB 
GMi=  40.25  dB 
PM  =1.6  deg 

Note  also  the  CLBW  is  t^proximately  3.5  rad/sec.  For  a  sample  rate  of  16  Hz,  the 
CLBW  would  need  to  be  about  25  rad/sec  before  it  would  violate  the  guidelines  discussed 
in  Section  (5.8). 


alpha  (deg)  de'^ae  (deg)  ^  _  Normal  g  (^) 
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6.3  F-I6  QDES  Design,  9.8  g  Response  Figures  (6.4)  and  (6.5)  show  the 
results  of  the  QDES  design  for  a  9.8  g  command.  This  value  was  chosen  as  it  appeared  to 
be  near  the  maximum  N^.  The  plots  of  input  sensitivity  and  CLBW  are  not  shown  for  the 
rest  of  the  F-16  examples  as  they  do  not  change  appreciably  as  the  g  level  is  increased. 
There  is  always  a  spike  of  approximately  magnitude  45  at  the  Phugoid  firequency,  and  the 
plot  is  constrained  below  25  after  that  spike.  The  CLBW  varies  slightly  between  1.5  and  4 
rad/sec  but  never  approaches  the  25  rad/sec  limit  previously  discussed.  Notice  that  the 
AOA  is  at  18  deg,  so  there  may  not  be  potential  for  a  higher  g  level.  Also,  although  the 
elevator  deflection  is  well  within  limits,  the  deflection  rate  needs  to  be  examined.  For  the 
F-16,  the  maximum  elevator  deflection  rate  is  60  deg/sec.  The  "burble"  in  the  Njj  response 
at  i  sec  is  caused  by  displacement  of  the  elevator  in  the  wrong  direction  (positive)  at  a 
rate  exceeding  100  deg/sec.  If  this  rate  were  required  to  achieve  the  desired  g  level  in  the 
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nose  up  direction,  this  value  would  be  of  concern.  However,  since  the  aircraft  is  going  the 
wrong  way  at  this  point,  the  designer  does  not  want  the  elevatOT  to  be  able  to  re^nd  that 
quickly  and  so  the  fact  that  it  cannot  is  not  of  concern.  In  the  time  period  between  1  and  2 
seconds  where  the  elevator  is  deflecting  to  cause  the  desired  level  of  performance,  the 
deflection  rate  is  approximately  29  deg/sec  which  is  well  within  limits. 
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6.3.3  F-16  QDES  Design,  10.0  and  10.2  g  Responses  This  section  compares  the 

responses  for  the  QDES  designs  for  10.0  and  10.2  g’s.  Comparing  figures  (6.6)  and  (6.7), 
nonce  how  QDES  has  eliminated  the  overshoot  from  the  10.0  g  response  in  order  to  keep 
the  AOA  within  limits  allowing  10.2  g's  to  be  achieved.  This  is  also  seen  when  coiiq)aring 
the  AOA  responses  in  figures  (6.8)  and  (6.9).  Note  that  for  10.0  g's,  the  AOA  limit  is  hit 
at  1.8  sec  corresponding  to  when  the  N2  response  hits  its  maximum  overshoot.  For  the 
10.2  g  case,  the  18  AOA  limit  occurs  at  about  3  sec  whe.i  the  response  is  within  the  3  % 
Ccmd  envelope. 
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Normal  g  (gs)  Normal  g  (gs) 
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Figure  6.8.  F-16  Elevator  and  AOA  10.0  g  Respcmses 
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6.3.4  F-16  QDES  Design,  10.4  g  Response  Hnally,  the  limits  of  performance  fw 

the  given  constraints  are  reached  at  10.4  g's.  Notice  how  the  N2  response  has  flattened 
out.  Also,  the  a  plot  stays  right  at  the  18  deg  limit  after  reaching  it.  The  maximum 
elevator  deflection  rate  (excluding  the  "burble")  is  32  deg/sec.  Figures  (6.10)  and  (6.11) 
show  these  results. 
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Figure  6. 11.  F-16  Elevator  and  AOA  10.4  g  Responses 


6.4  F-16  Short  Period  Approximation 

As  in  the  A-4D  example,  QDES  was  tried  on  the  Short  Period  approximation  for  the 

F-16.  Initially,  only  the  forward  velocity  (v)  state  was  removed  as  the  6  column  of  the  A 

matrix  did  not  become  all  zeros.  However,  this  left  a  low  ftequency  RHP  pole  at  .0015 

which  prevented  S  ftom  being  constrained  below  approximately  45.  Therefore,  both  the  v 

and  9  state  were  removed  leaving  the  following  poles: 

-20.00 

-4.350 

+1.905 

However,  this  did  not  completely  clear  up  the  sensitivity  constraint  problem  as  it  did  with 
the  A-4D.  QDES  was  not  able  to  return  a  solutiem  with  the  sensitivity  constrained  below 
20  which  resulted  in  the  following  margins  for  all  runs: 
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GMu=-0.4dB 
GMi=40.4dfi 
PM  =2.8  deg 

Constraining  Cemd  zero  was  also  tried  as  it  was  successful  in  the  A-4D  Short  Period 
case.  However,  for  the  F-16,  this  constraint  made  the  response  more  oscillatory  than 
without  it.  Therefore  it  was  not  used. 

6.4.1  F-16  Short  Period  QDES  Design,  1  g  Response  The  F-16  Short  Period 

QDES  design  response  to  a  1  g  command  is  shown  in  figures  (6.12),  (6.13),  and  (6.14). 
Notice  the  N2  response  is  smoother  than  the  full  model.  The  elevator  initially  deflects  the 
wrong  way,  resulting  in  some  nonminimum  phase  behavior.  The  sensitivity  plot  does  not 
have  a  low  frequency  spike  due  to  the  removal  of  the  Riugoid  mode;  however,  it  cannot 
be  constrained  below  tq)proximately  20.  The  closed  loop  bandwidth  is  iq)proximately  1.6 
rad/sec,  again  well  below  the  25  rad/sec  limit  discussed  in  Section  (6.3.1).  The  sensitivity 
and  CLBW  plots  are  not  shown  for  subsequent  Short  Period  runs  as  they  are  virtually 
identical  in  each  case. 


6-12 


alpha  (deg)  .<teltae  (deg)  Nomial  g  (gs) 


Figure  6.12.  F-16  Short  Period  QDES  Design,  1  g  Response 
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Figure  6.14.  F-16  Short  Period  1  g  Input  Sensitivity  and  CLBW 


<5.4.2  F-16  Short  Period  QDES  Design.  10.4  g  Response  Since  the  F-16  full 

model  maximum  was  10.4  g,  the  Short  Period  approximation  was  next  attempted  at 
this  level.  The  responses  are  shown  in  figures  (6.15)  and  (6.16).  Notice  that  due  to  the 
smoothness  of  the  response,  the  AOA  has  not  yet  hit  the  18  deg  limit,  even  though  the 
aircraft  is  delivering  a  10.4  g  response.  Also,  the  maximum  elevator  deflection  rate 
(excluding  the  nonminimum  phase  behavior)  is  only  19  deg/sec  contrasted  with  32  deg 
/sec  for  the  full  model.  The  fact  that  the  elevator  does  not  have  to  coiiq)ensate  for  a 
"burble"  at  approximately  1  sec,  as  did  the  full  model,  results  in  less  ctmtrol  usage. 
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alpha  (deg)  Nornial  g  (gs) 
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6 A. 3  F-16  Short  Period  QDES  Design,  10.6  g  Response 

The  limit  of  peifonnance  for  the  F-16  Short  Period  QDES  design  was  10.6  g‘s.  The 
responses  are  shown  in  figures  (6.17)  and  (6.18).  Notice  that  the  N2  response  hugs  the 
lower  3  %  e^md  boundary,  so  this  design  may  not  acnrally  be  delivering  a  higher  level  of  g 
than  the  full  model's  value  of  10.4.  From  an  operational  standpoint,  addition  of  0.2  g's  has 
little  significance.  However,  &om  a  strictly  analytical  point  of  view,  given  the  identical 
constraints  that  both  models  were  subjected  to,  diis  process  detemoined  that  the  Short 
Period  was  capable  of  a  higher  level  of  peifOTinance.  Of  significance  is  once  again  the 
lower  elevator  deflection  rate  (20  deg/sec)  required  by  the  Short  Period  design  to  achieve 
the  same  performance  level. 


Figure  6.17.  F-16  Short  Period  QDES  Design,  10.6  g  Response 
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Figure  6.18.  F-16  Short  Period  Elevator  and  AOA  10.6  g  Responses 


6.5  Summary 

This  chapter  explored  using  QDES  to  determine  the  limits  of  performance  for  the  F- 
16  as  was  done  in  Chapter  V  for  the  A-4D.  This  was  of  particular  interest  in  that  the  F-16 
is  a  much  more  difficult  aircraft  to  control  and  does  not  lend  itself  so  easily  to  a  classical 
design  as  did  the  A-4D.  The  ability  of  QDES  to  quickly  design  a  controller  for  varying 
performance  demands  proves  its  worth  in  a  situation  like  this.  However,  the  increased 
sampling  rate  required  for  the  F-16  did  increase  the  problem  size  and  slow  the  process 
down  considerably  compared  to  the  A-4D  example. 

The  F-16  full  model  responses  were  generally  srrxxither  than  the  A-4D.  This  may  be 
due  to  a  higher  number  of  FIR  taps  being  used  (45  vs.  35)  as  allowed  by  a  higher  number 
of  sample  times.  As  the  limiting  value  of  N2  was  reached,  it  appears  QDES  was  better 
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able  to  manipulate  the  time  response  of  the  F-16,  in  particular  the  overshoot,  thus  keeping 
AOA  from  hitting  the  limiting  value  of  18. 

The  Short  Period  approximation  responses  were  generally  smoother  than  their  full 
model  counterparts.  Both  achieved  a  slightly  hi^^  level  of  perftmnance,  due  to  the  less 
oscillatwy  nature  of  their  responses.  With  fewer  transients  in  Nx,  and  thus  the  a 
responses,  the  average  amplitude  of  the  AOA  response  was  higher  when  the  upper  limit  of 
its  oscillation  touched  up  against  the  18  deg  limit  This  allowed  slightly  nxne  elevatOT 
deflection  and  consequendy  slighdy  more  g. 

The  constraint  of  the  sensitivity  and  thus  control  of  the  stability  margins  was  not 
successful  for  the  F-16.  In  the  full  model  example,  the  spike  in  S  at  the  Phugoid  frequency 
remained  through  all  the  runs,  unlike  the  A-4D  example  where  it  disappeared  for  any 
response  over  1  g.  The  F-16  Short  Period,  althou^  lacking  the  low  frequency  spike, 
could  not  be  constrained  below  20.  This  resulted  in  unacceptably  low  stability  margins. 

The  author  is  aware  that  the  F-16  is  not  intended  to  be  flown  at  greater  than  9  g’s. 
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VII.  Step-By-Step  Procedure 


7.1  Introduction 

This  chapter  provides  a  detailed  step-by-step  procedure  for  using  QDES  to  evaluate 
aircraft  limits  of  performance  as  was  done  in  Chapters  V  and  VI.  The  intent  is  to  pipvide 
sufficient  detail  and  explanation  to  enable  the  non-controls  engineer  to  do  so.  Although 
the  explanation  will  use  the  evaluation  of  maximum  N2,  this  can  easily  be  adapted  to  any 
of  the  longitudinal  parameters.  It  could  also  be  adapted  to  the  lateral-directicmal  problem 
if  one  is  sufficiently  familiar  with  that  area.  To  avoid  needless  repetition,  this  chapter 
should  be  used  in  conjunction  with  this  diesis  and  its  appendices  as  well  as  the 
documentation  provided  with  the  QDES  program. 

7.2  QDES  Installation  and  Test 

Install  QDES  on  the  computer  system  per  the  instructions  in  the  QDES  user’s  guide 
supplied  with  the  tape.  It  will  be  necessary  to  copy  all  the  MATLAB  script  files  to  be 
used  into  the  user’s  home  directory  or  sub-directory,  wherever  the  work  is  to  be 
performed.  This  implies  that  the  user  must  have  MATLAB  available.  The  files  supplied 
with  the  tape  require  MATLAB  version  2.00  or  higher.  The  files  in  Appendix  C  were 
written  for  version  4.0.  The  author  found  it  useful  to  make  paper  copies  of  all  the 
README  files  supplied  with  the  individual  scripts  and  programs.  The  "Two  Mass" 
design  example  should  be  run  to  check  the  correct  operation  of  the  system.  Barratt  [2] 
offers  two  options  for  this.  The  use  of  '’qdes_set.m"  must  be  used  to  check  the  ability  of 
the  system  to  generate  the  impulse  response  matrices.  Primarily  tiiis  will  check  if  the  file 
’’mat2qdes"  is  installed  in  the  proper  path.  If  not,  a  message  will  be  generated.  TTiis  is 
explained  in  the  README  file  in  the  "sre"  sub-directory.  If  the  other  "Two  Mass"  option 
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is  chosen,  "mat2qdes"  will  not  be  called  and  QDES  will  use  impulse  response  matrices 
already  generated  and  supplied  with  the  program. 

Once  QDES'  correct  operation  has  been  verified,  the  following  MATLAB  files  must 
be  copied  into  the  directory  were  the  woric  will  be  performed.  Note  there  are  many  other 
script  files  supplied  with  the  QDES  program  tape;  however,  for  the  process  presented 
here,  these  are  the  only  ones  needed.  For  more  information,  see  [8].  These  fil^ are: 
fir2ss.m,  gethzw.m,  loadq.m,  midimpulse.m,  obc2ttt.m,  obc2ttt_fii.m,  qobc2comp.m.  It 
may  be  useful  to  make  paper  copies  of  these  files  for  easy  reference  should  a  problem 
occur. 

13  Model  Source 

The  A,  B,  C,  and  D  matrices  in  these  examples  are  derived  from  the  linearized 
equations  of  motion  for  the  aircraft  at  a  given  equilibrium  flis^ht  condition.  These  matrices 
can  be  formed  from  the  stability  derivatives.  Formulas  for  the  estimation  of  these 
derivatives  and  construction  of  the  matrices  can  be  found  in  Etkin  [9]  and  other  flight 
dynamics  textbooks.  It  is  critical  that  the  units  and  the  order  of  the  states  are  known. 

7.4  Using  the  "nzprpjn"  File 

As  briefly  described  in  Chapter  V,  a  MATLAB  file  similar  to  the  "nzprp.m"  shown 
in  Appendix  (C),  Section  (C.l)  should  be  written.  This  section  will  walk  through  this  file 
in  detail,  explaining  how  the  various  entries  are  arrived  at  and  the  rationale  behind  them 
based  on  the  author's  and  others  experiences  using  QDES. 

7.4.1  Entering  Plant  Data  The  plant  A,  B,  C,  and  D  matrices  are  entered.  To 
avoid  confusion,  it  is  recommended  that  the  order  of  the  states  match  those  used  in  this 
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thesis.  Next,  the  actuator  dynamics  are  augmented.  Appendix  (A),  Section  (A.  1.1)  and 
(A.  1.2)  show  how  the  matrices  change  as  the  state  space  representation  of  the  actuator  is 
augmented  to  the  system.  The  script  file  does  this  automatically.  Normally,  the  state- 
variable  model  matrices  are  constructed  such  that  N2  is  positive  nose  downward.  This  is 
easily  changed  to  the  positive  nose  up  convendon  used  in  this  diesis  1^  muldplying  the 
"nz"  Oast)  row  of  the  C  matrix  by  -1.  If  the  units  are  not  as  desired,  diey  can  be  ^isily 
converted.  Form  the  diagonal  T  matrix  to  scale  the  states.  The  number  on  the  diagonal  is 
what  must  be  multiplied  to  the  state  to  convert  fiom  the  present  units  to  the  desired  units. 
In  this  example,  the  original  states  are  in  non-dimensional  velocity  and  radians.  To 
convert  to  ft/sec  and  deg,  form  the  matrix: 


T  = 


'634 

0 

0 

0 

.  0 


0  0  0  0' 

57.3  0  0  0 

0  57.3  0  0 

0  0  57.3  0 

0  0  0  57.3. 


where  the  equilibrium  velocity  (V^)  is  634  ft/s^  and  1  rad  is  57.3  deg.  To  convert  the 
input  to  tleg,  form: 

S  =  [57.3] 


Finally,  to  convert  the  outputs  from  non-dimensional  velocity,  rad,  and  ft/sec,  form: 


634  0  0  0  0 

0  57.3  0  0  0 


W  = 


0 


0  57.3  0 


0 


0 

0 


0  0 
0  0 


57.3 

0 


0 

1 

32.2. 
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where  Ig  is  32.2  ft/sec.  Once  these  matrices  are  f(»med,  the  script  file  multiplies  them  to 
the  original  plant  matrices  as  shown  in  the  file  to  convert  them  to  the  new  units. 

As  mentioned  in  Chapter  V,  the  scale  facttrr  is  entered  to  tell  QDES  what  a  "unit  step" 
input  really  is. 

7.42  Discretizing  tite  Plant  OuqMer  V  provided  some  guidelines  to  selection  of 
the  sanofile  rate  (Srate  in  the  file)  to  be  entered.  The  author  recommends  using  the 
discretization  rate  based  on  five  times  the  frequency  of  the  fastest  system  pole  (usually  the 
actuator).  Use  of  a  lower  rate  risks  not  capturing  all  the  plant  dynamics.  In  some  cases, 
use  of  a  higher  rate  may  result  in  a  smoother  response  as  in  the  case  of  the  A-4D  Short 
Period  example,  so  there  is  room  for  expeiimentaticm.  The  author  found  that  increasing 
sampling  rate  tended  to  raise  the  minimum  value  below  which  the  sensitivity  could  be 
constrained.  Also,  increasing  the  sampling  rate  greatly  increases  the  size  of  the  problem 
with  its  attendant  slowing  of  the  process  as  was  experienced  with  the  F- 16  in  Ch^ter  VI. 

7.4.3  Forming  the  P-K  Style  Matrices  The  formation  of  the  nine  matrices  needed 
by  "obc2ttt.m"  is  straightfinward  as  developed  in  Chapter  V.  These  can  be  used  exactly 
as  depicted  in  the  file  unless  it  is  desired  to  change  the  number  or  order  of  the  various 
inputs  or  outputs  (w,  z,  y,  u).  If  so,  a  diagram  such  as  figure  (5.1)  must  be  drawn  and  the 
matrices  constructed. 

7.4.4  Forming  an  Observer  Based  Controller  As  explained  in  Chapter  V,  QDES 
must  be  given  a  stable  nominal  plant  with  a  nominal  controller.  In  this  process,  this  is 
accomplished  by  forming  an  Observer  Based  Controller  using  the  MATLAB  "dlqr"  and 
"dlqe"  cennmands  to  generate  the  state  feedback  gain  matrix  (K)  and  the  Kalman  filter  gain 
matrix  (L).  It  is  recommended  that  the  wdghts  listed  (Chapter  V  and  as  entered  in  the 
example  file  be  used.  The  author  experimented  with  other  weights  and  with  using  the 
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MATLAfi  "place"  command  as  done  in  the  “Two  Mass"  exan^le.  Both  cases  produced 
worse  results.  In  theory,  as  long  as  the  nominal  system  is  stable,  it  shouldn’t  matter  to 
QDES  what  it  is.  The  resulting  solutitxi  should  be  the  same.  This  was  certainly  not  the 
case. 


7.45  Calling  "obc2ttt.m"  Two  other  pieces  of  information  need  to  be.  entered 
before  calling  "obc2ttt.m"  to  form  the  impulse  response  matrices.  The  first  is  the  source 
file  name.  This  name  is  used  to  tie  the  QDES  source  file,  the  impulse  response  matrices, 
and  the  output  files  together.  Whatever  is  entered  here  will  be  tacked  on  to  several  files 
associated  with  the  problem  being  addressed  by  the  QDES  source  file.  The  other  piece  of 
information  is  the  number  of  samples  (n_sample).  This,  along  with  the  sample  rate, 
determines  over  what  time  span  QDES  will  search  for  a  solution.  It  is  the  sample  rate  in 
Hz  times  the  time  over  which  one  wishes  QDES  to  work.  As  noted  in  Chapter  VI,  the 
number  of  samples  relates  directly  to  the  si^  of  the  problem  in  terms  of  computer 
memory.  Therefore,  this  number  should  be  kept  as  small  as  possible.  Widi  a  fixed  sample 
rate,  this  means  the  time  span  should  be  kept  as  small  as  possible.  As  discussed  in  Section 
(5.4),  for  the  N2  problem  there  is  litde  point  in  running  the  problem  for  a  long  period  of 
time.  Three  to  five  seconds  should  be  reasonable  in  most  cases,  unless  the  aircraft  has  an 
unusually  slow  response.  Assuming  there  are  no  errors,  "obc2ttt.m"  writes  the  following 
files  to  the  current  directory;  A4D7NZ_tl,  A4D7NZ_t2,  A4D7NZ_t3,  and 
A4D7NZ_ttt.mat  Note  these  use  the  source  fife  name  fiom  the  exanq)le  script.  The  first 
three  are  the  impulse  response  matrices,  and  the  last  ctmtains  the  plant  information  and 
nominal  controller  for  use  by  QDES  and  for  forming  the  closed  loop  plant  once  QDES 
arrives  at  a  solution.  The  "nzpip.m",  as  written,  also  displays  the  following  message  on 
the  screen; 

When  writing  your  Qdes  source  file,  be  sure  that 


7-5 


njsanple-56 

so  that  it  inatdies  with  the  iiqpulae  ze^xnse  matrices  3/ou  just  created 

For  this  set  up  your  sanple  rate  was:  16  Hz  Sanple  tine  was:  0.0625  sec 
Also,  be  sure  that  your  "nzpst  jn"  file  uses  the  Scale  Factor 
S£nz-  1.0 

7  J  VVrin'ng  the  QDES  Source  File 

This  sectim  walks  through  the  writing  of  the  QDES  source  file.  The  purpose  and 
use  of  this  file  was  discussed  in  Chapter  V.  This  sectitm  will  detail  the  mechanics  of 
putting  it  together.  Appendix  (C),  Section  (C.2)  shows  the  source  file  used  for  the  A-4D 
example  and  will  be  referenced.  This  file,  written  in  "C",  must  be  given  the  same  name  as 
that  entered  under  "source  in  the  "nzprp.m"  script  file. 

7.5./  The  Preamble  In  the  preamble,  the  exogenous  inputs  (w),  regulated 
outputs  (z),  and  sensed  outputs  (y)  are  defined  and  ordered  or  numbered.  This  order  must 
match  that  used  when  constructing  the  matnces  in  "nzprp.m”.  When  writing  the 
constraints,  these  signals  may  be  referenced  either  by  their  name  as  defined  or  their 
number.  In  this  example,  numbers  were  used  exclusively. 

75.2  The  Declarations  Section  In  the  declarations,  the  number  of  each  of  these 
signals  listed  in  the  preamble  must  be  listed,  as  well  as  the  number  of  acmators.  The 
number  of  samples  from  the  "nzprp.m"  file  is  entered  as  well  as  the  number  of  FIR  taps  to 
be  used.  As  mentioned  in  Chapter  V,  the  number  of  taps  is  constrained.  The  author  found 
a  good  rule  of  thumb  was  to  make  the  number  of  taps  five  less  than  the  number  of  san^les 
unless  this  would  cause  the  number  of  taps  times  number  of  actuators  times  number  of 
sensed  outputs  to  exceed  200,  or  if  it  made  the  problem  too  large  in  terms  of  computer 
memory.  The  rationale  for  using  as  many  taps  as  possible  is  that  the  goal  is  to  see  if  a 
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controller  can  be  built  to  meet  the  constraints.  If  more  taps  are  used,  QDES  has  more 
information  with  which  to  work  and  can  better  answer  that  question.  In  nK>st  cases,  using 
a  hi^er  number  of  taps  also  resulted  in  a  smoother  response.  Hnally,  the  inpilse 
response  matrices  must  be  defined  by  name.  The  example  source  file  can  be  used  as  a 
model  and  should  be  followed  exacdy  to  ensure  proper  syntax. 


where:  DUGM  =  Desired  Upper  Gain  Margin  (dB) 

DLGM  =  Desired  Lower  Gain  Margin  (dB) 

DPM  =  Desired  Phase  Margin  (deg) 

Note  that  only  one  value  is  entered  (the  lowest)  idthough  there  are  three  separate 
constraints.  Thus,  whatever  constraint  gives  the  lowest  value  drives  the  ||Sy^.  If 


achieved,  the  other  two  margins  will  be  bett^  dian  desired.  On  the  other  hand,  if  QDES 
misses  on  constraining  |sy„,  this  may  not  effect  two  of  the  margins,  depoiding  on  how 
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much  it  missed  by.  When  using  this  constraint  function,  QDES  warns  the  user  that  a  0.34 
dB  eiTor  may  result.  This  equates  to  a  factor  of  1.0399.  Thus,  if  a  certain  value  of  is 

absolutely  required,  this  value  should  be  multiplied  by  1.0399  and  the  difference  between 
this  and  the  desired  value  should  be  subtracted  from  the  desired  value  and  entered  in  the 
constraint.  This  assumes  QDES  will  error  on  the  high  side  which  would  give  lower 
margins  than  desired  if  not  compensated  for.  Ccmstraining  the  maximum  elevator 
deflection  and  angle  of  attack  are  done  via  the  ''u_step[  ][  ]  (t)"  function.  Again,  the 
output  is  listed  first,  then  the  input.  Note  that  the  time  (t)  is  discrete  time  steps,  not  actual 
clock  time.  The  elevator  constraint  in  the  source  file  says  that  for  a  unit  step  input  to  [1] 
(Njcntd)  output  [4]  (S^)  should  stay  within  the  values  listed.  The  same  idea  holds  for 
the  angle  of  attack  constraints.  Note  that  needs  to  be  taken  into  account  when 
setting  the  numbers  in  the  constraint  as  noted  in  Section  (S.S.4).  The  overshoot  and 
undershoot  constraints  are  fairly  straightforward.  To  draw  the  remaining  time  dmnain 
constraint  curves  and  lines  shown  in  response  plots  in  Qiapters  V  and  VI,  various 
exponential  relationships  are  required.  To  accomplish  this,  a  MATLAB  script  file 
"visnz.m"  was  written.  This  file  is  found  in  Appendix  (C),  Section  (C.4).  By  following 
the  comments,  one  can  enter  the  desired  time  domain  constraints  and  the  program  will 
figure  all  the  exponential  relationships  and  draw  a  sample  envelope  for  examination.  By 
varying  the  inputs,  one  can  tailor  the  envelope  to  virtually  any  shape  desired.  Most 
importantly,  the  script  file  also  writes  out  the  constraints  as  they  should  appear  in  the 
QDES  source  file.  An  example  of  the  output  from  ’’visnz.m"  for  the  A-4D  runs  is  shown. 

This  is  how  your  Qdes  Source  File  constraint  statenents  should 
look  to  draw  the  equivalent  constraint  envelope  for  Qdes  and  constrain 
the  vector  margins  of  your  selected  SISO  transfer  function. 

***************  jxxt  NCPMAL  AXEIERATICN  (Nz)  ********** 


/*  Put  Constraints  On  "  "  Vector  Margins  */ 

for  %(-0.0  to  %pi  st^  0.02  do  { 
nBg_H[out]  [in]  (l,w) <-1.663; 

) 

/*  Put  Constraints  On  Overshoot  */ 
overshoot [out] [in]<=  0.1; 

/*  Put  Ccsistrains  On  Undershoot  */ 
undershoot [out] [in]<-0 . 05; 

/*  Put  Constraints  On  Settling  Tine  */ 
for  t-l6  to  20  do 

l-0.8392^<-u_step[out]  [in]  (t)<-l+0.8392'^; 

/*  Request  Zero  Ccmmand  Error*/ 

Re_H[out]  [in]  (1,0)  =  1; 

MO/OR 

/*  Put  Constraints  On  Ccmtand  Error  */ 
for  t=20  to  n_3anple-l  do 

0.97+0* (t-20)<-u_step [out]  [in]  (t)<-1.03-0*(t-20) ; 

/*  Put  Constraint  On  Rise  Tine  */ 

for  t=12  to  n_sairple-l  do 
0.8<=u_st^[out]  [in]  (t) ; 

/*  Put  In  l^jper  Constraint  Curve  (if  desired)  */ 
for  t=0  to  4.8  do 

u_step[out]  [in]  (t)<=(0.1596*t)  ''3+0.2; 
for  t=4.8  to  12  do 

u_st^[out]  [in]  (t)<=0.45-(0. 7663-0. 1064* (t-4.8))"3+0. 65; 
/*  Put  In  Lower  Constraint  Curve  (if  desired)  */ 
for  t-4  to  7.2  do 

(0.2303*  (t-4) )  ''3<=u_step[out]  [in]  (t) ; 
for  t=7.2  to  12  do 

0.4-(0. 7368-0. 1535*(t-7.2))''3+0.4<=u_step[out]  [in]  (t); 
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In  a  windows  environment,  it  is  then  a  simple  matter  of  cutting  and  pasting  the  ly^propriate 
constraint  to  the  QDES  source  file.  The  only  other  thing  that  must  be  done  is  to  insert  the 
name  or  number  of  the  output  and  input  of  interest  in  the  brackets  where  it  currently  says 
[out]  [in]  and  double  check  the  syntax  for  proper  brackets  as  shown  in  the  exan^le  source 
file. 


7.6  Running  QDES 


To  run  QDES,  type  "qdes  A4D7NZ"  on  the  command  line  and  hit  ittum.  If  QDES 
finds  any  syntax  or  other  errors  it  wiU  stop  and  provide  a  message.  These  are  fairly  clear, 
identifying  what  line  of  the  source  file  contains  the  error.  Other  times  it  may  identify  a 
mismatch  between  the  number  of  samples  entered  in  the  source  file  and  the  number 
actually  in  the  impulse  response  matrices.  Again,  these  are  fairly  obvious.  If  there  are  no 
errors,  QDES  returns  the  following: 

qdes  version  1.5  infonnal  distribution 
cxjpyright  1987  Stanford  University 
got  all  (xefficients 

Warning:  Feferenoe  to  Mag_H  may  incur  +/-0.34  dB  error. 

329  control  constraints 
0  (±ijective  terms 
generating  feasibility  program 


Calls  to  LSCFm 


print  level  »•  10 
feasibility  tolerance  *  le-06 
crash  tolerance  =  0.0001 
infinite  bound  size  =  le+10 
print  level  =  0 
solving  feasibility  program 
problan  type  FP 

If  QDES  finds  a  solution,  the  following  will  be  letumed: 

found  feasible  point  in  62  iterations 
writing  F162NZ.q  F162NZ.out 
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If  no  solution  is  found  QDES  returas  the  following: 

no  feasible  point  found  in  131  iterations 

svm  of  infeasibilities  >  0.0459192 

writing  F162NZ.last.q  F162NZ.out 

When  the  sum  of  the  infeasibilities  is  quite  small  (less  than  0.1),  adding  a  few  more  taps  (if 
possible)  may  result  in  a  solution.  The  ".out"  file  may  be  examined  to  see  what  constraints 
QDES  could  not  meet.  For  multiple  runs,  if  the  sample  rate,  number  of  samples,  or  scale 
factor  have  not  changed,  "nzprp.m"  does  not  have  to  be  ran  again.  Constraints  in  QDES 
can  be  changed  and  more  runs  attempted.  If  the  sample  rate  changes  however,  "visnzjn" 
will  have  to  be  ran  with  the  new  value  as  all  the  calculations  for  the  constraint  curves  are 
based  on  this.  When  Erst  starting  a  QDES  problem,  it  is  reccHomended  that  the 
constraints  be  added  in  progressively.  Otherwise,  if  no  solution  is  found,  it  may  be 
difficult  to  ascertain  which  constraint  is  causing  the  difficulty.  The  ej^j  constraint  should 
be  entered  first  to  get  QDES  to  form  a  step  re^nse.  The  other  time  constraints  can  then 
be  added  starting  from  the  end  of  the  response  (T^i)  and  woridng  forward  (Tj).  The 
margin  constraint  should  be  entered  after  all  the  time  constraints  have  been  explored  as  it 
tends  to  lengthen  the  ran  times.  Elevator  and  AOA  constraints  should  be  added  only 
when  the  limits  of  performance  are  to  be  explored,  as  they  greatly  increase  the  problem 
size  and  slow  the  confutations.  Use  of  the  zero  command  error  constraint  should  be 
reserved  until  all  other  objectives  are  met.  It  may  improve  the  response  as  it  did  for  the  A- 
4D  Short  Period,  but  more  often  than  not,  it  may  result  in  no  solution  being  returned. 

7.7  Outputting  the  Results 

The  MATLAB  script  file  "nzpst.m"  in  Appendix  (C),  Section  (C.3),  takes  the  QDES 
solution,  forms  the  closed  loop  system,  and  outputs  the  information  for  each  ran  as  seen  in 
(Chapters  V  and  VI.  The  QDES  source  file  name  must  be  entered  to  enable 
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"qobc2coinp.in"  to  form  the  QDES  controller  and  "gethzwm"  to  form  the  closed  loop 
system.  The  sample  rate,  run  time  (sec),  and  scale  factor  are  entered  and,  of  course,  must 
match  those  used  to  set  up  the  problem.  The  data  horn  "visnzm"  should  also  be  entered 
to  enable  the  file  to  draw  the  same  constraint  envelope  as  that  imposed  on  QDES.  In  the 
section  of  the  file  where  the  N2  time  response  is  plotted,  each  plot  command  for  the 
various  time  constraint  curves  is  on  a  separate  line  and  labeled.  These  may  be  individuaUy 
commented  out  as  desired  (using  %)  to  plot  c«ly  those  constraints  which  were  actually 
imposed  on  QDES.  To  evaluate  elevator  deflection  rates,  another  MATLAB  script  file 
"surf.m"  was  written  and  is  found  in  Appendix  (C),  Section  (C.5).  The  file  contains  simple 
instructions  for  its  use.  When  used,  it  returns  the  surface  deflection  rate  over  the  time 
frame  designated. 

7.8  Summary 

This  chapter  developed  a  simplified,  step-by-step  approach  to  using  QDES  to 
evaluate  aircraft  N2  limits  of  performance  as  accomplished  in  Chapters  V  and  VI. 
Although  not  completely  self-contained,  when  used  with  the  body  of  this  thesis  and  its 
appendices,  it  provides  most  of  the  nuts-and-bolts  information  needed  by  a  non-controls 
type  engineer  to  perform  this  type  of  analysis. 
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V7//.  Conclusions  and  Recommendations 


8.1  Introduction 

This  chapter  presents  an  evaluation  of  using  QDES  to  explore  the  limits  of 
performance  in  aircraft  design  as  was  acctxnplished  in  this  thesis.  It  also  offers  possible 
explanations  as  to  why  some  aspects  of  the  process  did  not  work  as  well  as  anticipated  as 
well  as  other  approaches  tried.  Hnally,  it  offers  a  suggestion  for  future  research  in  this 
area. 


8.2  Usefulness  of  QDES 

QDES'  method  of  operation  is  quite  ingenious.  The  ability  to  enter  a  number  of 
constraints  and  then  let  the  computer  see  if  it  is  possible  to  build  a  controller  to  meet  them 
lends  itself  perfectly  to  the  limits  of  performance  problem.  Evaluation  of  QDES’ 
usefulness  for  this  problem  can  be  divided  into  two  areas,  time  response  and  frequency 
response. 

8.2.1  Time  Response  Examining  the  time  response  portion  of  the  problem,  that 
is,  the  aircraft  N2,  a,  and  6^  responses  to  a  step  input,  QDES  did  quite  well.  Although  the 
Nz  responses  were  not  always  picture  perfect,  they  remained  within  the  constraint 
envelope.  Furthermore,  QDES  was  easily  able  to  manipulate  them  by  reducing  the 
overshoot  or  oscillation  of  the  response  to  accommodate  higher  Nzcmd  and  still  remain 
within  the  AOA  constraint.  This  capability  allowed  quick  determination  of  the  maximum 
achievable  performance  of  the  aircraft  at  the  specified  flight  condition.  The  results  from 
this  type  evaluation  will  indicate  whether  or  not  the  aircraft  can  be  controlled  up  to  the 
limits  of  performance;  however,  the  actual  values  of  N^,  a,  or  bg  should  be  viewed  with 
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caution.  This  is  because,  even  at  moderate  g  levels  over  a  short  time  period,  the  angles 
which  were  assumed  small  in  the  equations  of  motion  linearization  process  are  no  longer 
small.  This  is  a  limitation  of  the  modeling  of  the  problem  and  occurs  whenever  the 
linearized  equations  are  used  regardless  of  the  controller  design  methodology. 

8.2.2  Frequency  Response  The  success  of  QDES  in"ranstraining  the  input 
sensitivity  and  thus  evaluating  the  achievable  vector  stability  margins  was  mixed.  For  the 
A-4D  full  model  example,  the  sensitivity  had  a  high  amplitude  spike  at  the  Phugoid 
frequency  that  could  not  be  constrained  at  1  g  even  though  QDES  said  it  had  found  a 
feasible  solution.  The  same  phenomena  occurred  in  the  F- 16  full  nxxlel  example.  The 
most  logical  explanation  for  this  is  that  QDES  did  not  get  enough  plant  data  to  see  and 
control  the  Hiugoid  mode.  Although  the  sampling  rate  used  was  more  than  sufficient  to 
capture  the  low  frequency  Phugoid  (it  was  based  on  the  much  faster  actuator  pole),  the 
sampling  time  was  not  a  sufficient  fraction  of  the  relatively  long  (100  sec)  Phugoid  period. 
Thus,  the  impulse  response  matrices,  from  which  QDES  gets  its  infonnation  ,  did  not 
contain  enough  samples  for  QDES  to  see  the  Phugoid.  However,  this  does  not  explain 
why  the  sensitivity  was  constrained  above  1  g  for  the  A-4D  example  but  not  for  the  F-16. 
The  author  feels  this  has  to  do  with  the  information  given  to  QDES.  QDES  does  not  use 
the  uncompensated  plant  for  its  design  as  is  done  classically.  Rather,  it  takes  the  complete 
closed  loop  system  with  a  nominal  controller  already  in  place.  It  is  felt  that  differences  in 
this  nominal  plant  can  cause  different  results  from  QDES  in  the  area  of  the  input 
sensitivity.  For  example.  Appendix  (A),  Section  (A.2)  contains  the  state  feedback  gain 
matrix  (K),  the  Kalman  filter  gain  matrix  (L),  and  the  discrete  poles  of  the  ncmoinal  system 
for  each  A-4D  full  model  run.  Notice  that  although  K  remains  the  same  for  all  runs,  L  and 
the  discrete  poles  do  change.  The  most  significant  changes  in  the  discrete  poles  occur 
between  1  and  3  g's  when  the  sensitivity  goes  from  being  unconstrained  to  constrained. 
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From  3  g's  on  up,  the  values  are  similar,  as  are  the  fiequency  response  results.  The  F-16 
full  model  Appendix  (B),  Section  (B.2)  did  not  exhibit  exactly  the  same  behavitn*.  Only 
the  first  discrete  pole  changed  significantly  from  1  to  9.8  g's  and  then  remained  essoitially 
unchanged  for  the  rest  of  the  runs.  All  the  other  discrete  poles  remained  virtually 
unchanged  through  all  the  runs.  Interestingly,  the  low  fiequency  I%ugoid  poles  are  not 
the  ones  changed  when  going  from  1  g  to  higher  g  levels  in  either  the  A-4D  or  F-16  case. 
One  would  expect  that  for  a  change  in  the  ncaninal  system  poles  to  improve  the  QDES 
results,  this  change  would  be  a  movement  of  the  Phugoid  poles  to  a  higher  frequency. 
Turning  to  the  Short  Period  examples,  the  A-4D  Short  Period  performed  well  in  the 
fiequency  response  area  with  no  unconstrained  low  fiequency  spikes.  This  is  what  one 
would  expect,  having  removed  the  low  fiequency  Phugoid  mode.  However,  the  F-16 
Short  Period  did  not  perform  as  well.  Although  no  low  fiequency  spike  was  observed,  the 
sensitivity  across  the  entire  fiequency  range  would  not  be  constrained  below 
approximately  20  regardless  of  the  g  level.  Examination  of  the  change  in  the  discrete 
poles  of  the  nominal  system  from  1  to  higher  g  levels  show  the  same  behavior  as  the  F-16 
full  model.  The  first  pole  changes  significantly  fiom  1  g  to  the  higher  levels,  but  is 
unchanged  after  that.  The  remainder  of  the  poles  are  unchanged  throughout  all  the  runs. 
None  of  these  discrete  poles  are  low  frequency  Phugoid  type. 

8.3  Other  Approaches 

In  an  attempt  to  see  if  some  of  these  fiequency  problems  could  be  side-stepped, 
several  different  methods  were  tried  with  mixed  results. 

8.3.1  Longer  Run  Times  In  an  attempt  to  capnue  the  elusive  Phugoid  mode, 
some  examples  with  longer  run  times  were  tried  .  However,  to  get  a  solution,  the 
sampling  rate  had  to  be  reduced  to  well  below  the  minimum  required  to  capture  the  fastest 
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pole.  Even  then,  in  only  one  case  with  the  A-4D  sampled  at  3  Hz  for  20  seconds  was  the 
sensitivity  constrained.  This  combination  of  san^Ie  rate  and  nm  time  is  incompatiUe  with 
the  limits  of  performance  procedure  developed  in  this  thesis.  Aside  from  possibly  not 
capturing  all  the  aircraft  plant  dynamics,  the  use  of  a  long  run  time  would  result  in  no 
solutions  because  a  and  5^  would  increase  rapidly  as  forward  velocity  decreased  in  the 
pitch  up.  The  limits  would  soon  be  exceeded  with  no  QDES  solutitm  returned. 

8.3.2  Different  Nominal  Plant  Since  the  A-4D  full  model  responded  favorably  to 
a  change  in  the  iKxninal  plant  offered  to  QDES,  an  attempt  was  made  to  use  the 
MATLAB  "place"  command  to  generate  the  K  and  L  matrices  as  was  done  in  the  "Two 
Mass"  example  [2].  The  intent  was  to  move  the  Phugoid  poles  so  as  to  increase  their 
frequency.  Several  different  placements  were  tried  with  results  that  were  considerably 
worse  in  terms  of  the  frequency  respcmse  as  compared  to  the  weighing  used  in  the 
examples  shown  in  this  thesis. 

8.3.3  Output  Margins  An  alternative  to  constraining  and  observing  the  input 
sensitivity  is  to  do  so  at  the  output  of  the  plant.  This  requires  constraining  the  individual 
SISO  transfer  functions 

Nz 

Nzcmd 

a 

^zcmd 

q 

^zcmd 

This  was  attempted  on  the  F-16  fidl  model  and  Short  Period  approximations.  The  F-16 
full  model  output  sensitivities  behaved  in  the  same  manner  as  the  input  sensitivity.  The 
sensitivity  was  flat  at  unity,  while  the  a  and  q  responses  had  magnitude  spikes  at  the 
Phugoid  frequency  which  would  not  be  constrained  at  either  1  g  or  higher  levels.  The 
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Short  Period  output  sensitivity  were  constrained  below  2  in  all  cases.  Thus,  in  the  F- 16 
Shtfft  Period  case,  observing  the  margins  at  the  output  resulted  in  a  markedly  diffoTOt  and 
better  result  than  observing  the  input  mar^s. 

S.4  Conclusion 

The  limits  of  performance  methodology  presented  in  thesis  could  be  used  with  some 
cautions.  The  N2  performance  results,  subject  to  the  a  and  5^  constraints,  ate  reliable 
with  the  caveats  offered  in  Section  (8.2.1)  concerning  the  linearization  process.  However, 
the  results  from  the  vector  margin  evaluation  are  too  inconsistent.  Other  design  processes 
have  produced  controllers  for  the  F-16  with  phase  margins  near  60  deg.  An  area  for 
future  research  might  be  investigation  of  the  relationship  between  the  nominal  controller 
offered  to  QDES  and  the  frequency  results  returned.  To  be  truly  useful  for  aircraft  limits 
of  performance  evaluation,  QDES  needs  to  be  irxxliEed.  More  speed,  plus  continuous 
time  domain  and  MIMO  independent  stability  margins  capability  would  go  a  long  way  in 
overcoming  the  difficulties  encountered  in  this  thesis. 
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Appendix  A.  A-4D  State  Space  Models  and  Design  Information 


This  sq)pendix  contains  the  continuous  state  space  model  for  the  A-4D  at  flight 
condition  5.  This  model  was  used  for  the  classical  design  of  Chapter  IV  and  is  the  starting 
point  for  the  QDES  design  of  Chapter  V.  It  also  contains  the  state  feedback  gain  matrix 
(K),  the  Kalman  filter  gain  matrix  (L),  and  the  discrete  ptdes  of  the  nominal  closed  loq> 
system  provided  to  QDES  for  each  design  run.  The  matrices  and  vectws  are  as 
represented  in  the  state  variable  differential  equations 

x  =  Ax  +  Bu 
y  =  Cx+Du 

A .  1  A-4D  Longitudinal  Model  at  Flight  Condition  5 

A. 1.1  Original  Plant  Model  The  state,  input,  and  output  vectors  are  as  follows: 


’  V  ' 

v' 

a 

x  = 

a 

6 

u=[Se]  y= 

e 

q 

.q. 

Nz. 

where 

V  is  forward  velocity  {normalized  to  V^y 
a  is  angle  of  attack  {rad). 

6  is  pitch  angle  {rad). 
q  is  pitch  rate  {radtsec). 

N2  is  normal  acceleration  at  the  aircraft  center  of  gravity  ifttsec^). 
hf.  is  elevator  deflection  {rad). 
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The  matrices  associated  with  this  model  are  as  follows: 


-0.01266 

-0.00588 

-0.05075 

0  " 

0 

-0.10100 

-0.81670 
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0.99840 

-0.08963 

A  = 
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B= 

0 
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1 
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0 

0 
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1 

0 
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0 

0 

0 

0 

1 

0 
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-517.788 

0 

-1.02708, 

-56.8254 

A.I.2 

dynamics 


A-4D  Plant  with  Actuator  Dynamics 

6e  10 

Secmd  s+10 


Inserting  the  elevator  actuator 


represented  as 

A  =  [-10]  B  =  [10]  C=[l]  D*[0] 

and  changing  N2  to  positive  nose  up  yields 
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The  state,  input,  and  output  vectors  are  now 
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V 

V  ’ 

a 

a 

6 

u=[8ecmdl  y= 

e 

q 

q 

5c. 

.Nz. 

where 

V  is  fcHward  velocity  {normalized  to  V^). 
a  is  angle  of  atuuk  (rad). 

6  is  pitch  angle  (rad). 
q  is  pitch  rate  (radfsec). 

N2  is  normal  acceleration  at  the  aircraft  center  of  gravity  (fUsec^). 

is  elevator  deflection  (rad). 

5ecind  is  elevator  deflection  commanded  (rad). 

A.1.3  A-4D  Model  in  Desired  Units  Converting  forwaid  velocity  to  ft/sec,  angles 
and  angular  rates  to  degrees,  and  N2  to  g's  yields 
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0 

with  the  same  state,  input,  and  output  vectors  and  the  following  units 
V  is  forward  velocity  (ft! sec). 
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a  is  angle  of  attack  {deg). 

0  is  pitch  angle  {deg). 
q  is  pitch  rate  {deg! sec). 

Nz  is  normal  acceleration  at  the  aircraft  center  of  gravity  {g’s). 
Se  is  elevator  deflection  {deg). 

S^cmd  cicvator  deflection  commanded  {deg). 


AJ  A-4D  QDES  Design  Information 

This  section  ccmtains  the  feedback  gain  matrices,  Kalman  filer  gain  matrices,  and 
nominal  system  discrete  poles  for  the  various  QDES  design  runs  on  the  A-4D. 


A.2.1  A-4D  1  g  Command 


K  =  [  0.01311  0.19535 

-.01405 

-0.03883 

0.14333 

2,80394 
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2.63160 

0 

L  = 
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2.84941 
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Discrete  poles: 


-0.45308 

0. 74824  ±  0.34392i 

0.30119 

0.29055 

0.99892  ±0.007471 
0.98908  ±0.009011 
0.28650 


A.2.2  A-4D  3  g  Command 

K  =  [  0.01311  0.19535  -.01405  -0.03883  0.08394] 


A-4 
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0.07557 

5.78362 

0 
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A-4D  4J  g  Command 

K=[  0.01311 

0.19535 

-.01405 

-0.03883 

0.08394] 

-0.10271 

6.88706 

0 

-0.01630 

6.32552 
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Discrete  poles: 
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0.29055 

0.98908  ±0.00901i 
0.28650 

A.2.4  A-4D  4.6  g  Command  The  values  for  the  4.6  g  command  are  virtually 
identical  to  the  4.5  g  values. 
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A.3  A-4D  Short  Period  Approximation  Flight  Condition  5 


The  state,  input,  and  output  vectws  are  as  follows: 


x  = 

a" 
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II 

CO 

>1.1—1 

II 
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q 

8c. 

.Nz. 

where 

a  is  angle  of  attack  (deg). 
q  is  pitch  rate  (degisec). 

N2  is  nonnal  acceleration  at  the  aircraft  center  of  gravity  (g's). 
is  elevator  deflection  (deg). 

The  matrices  associated  with  this  model  are  as  follows: 
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A  .4  A-4D  Short  Period  Approximation  QDES  Design  / nformation 

This  section  contains  the  feedback  gain  matrices,  Kalman  filer  gain  matrices,  and 
nominal  system  discrete  poles  for  the  various  QDES  design  runs  on  the  A-4D  Short 
Period  approximation. 

A. 4.1  A-4D  Short  Period  1  g  Command 

K  =  [  0.06590  -0.03256  0.07101] 
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Discrete  poles: 
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A.4.2  A-4D  Short  Period  3  g  Command 

K  =  [  0.06590  -0.03256 


L  = 
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Discrete  poles; 
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A.4.3  A-4D  Short  Period  4.6  g  Command 

K  =  [  0.06590  -0.03256  0.07101] 


L  = 
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Discrete  poles: 


A-7 
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AAA  AAD  Short  Period  4.7  g  Command 

K  =  [  0.06590  -0.03256  0.07101] 
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Discrete  Poles: 
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Appendix  B.  F-16  State  Space  Models  and  Design  Information 


This  appendix  contains  the  continuous  state  space  model  for  the  F-16  at  sea  level 
and  Mach  0.6  as  used  in  Chapter  VI.  It  also  contains  the  state  feedback  gain  matrix  (K), 
the  Kalman  filer  gain  matrix  (L),  and  the  discrete  poles  of  the  mminal  ck>sed  loop  system 
provided  to  QDES  for  each  design  run.  The  state,  input,  and  output  vectOTS  are  the  same 
as  used  for  the  A-4D  in  Appendix  (A).  Initial  units  are  slightly  different  as  noted. 


B.l  F-16  Full  Longitudinal  Model 


The  matrices  associated  with  the  original  plant  without  actuator  dynamics  are: 
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-0.00214' 

A  = 

-0.00008 

0 

-1.49100 

0 

-0.00130 

0 

O.r  >600 

1 

B= 

0.18800 

0 

-0.00036 

9.75300 

0.00029 

-0.96000 

,  19.0400 

1 

0 

0 

0 

- 

■  0  ■ 

0 

1 

0 

0 

0 

C  = 

0 

0 

1 

0 

D= 

0 

0 

0 

0 

1 

0 

-0.00166 

-31.0239 

-0.027049  -0.08323  _ 

3.91 18_ 

with  units: 


v  (ft! sec) 
a  {rad) 

0  {rad) 
q  {rad! sec) 
Nz  {g's) 

5c  {rad) 
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Adding  the  actuator  dynamics: 

8e  -  20 
^ecmd  s+20 

as  in  Appendix  (A),  converting  to  the  same  units  as  used  for  the  A-4D,  and  changing  Nj 
to  positive  nose  up  yields 


■-0.01485 

0.65239 

-0.56195 

-0.31308 

-0.00004' 

01 

-0.00458 

-1.49100 

-0.00130 

0.99600 

0.18800 

0 

A  = 

0 

0 

0 

1 

0 

B= 

0 

-0.02063 

9.75300 

0.00029 

-0.96000 

19.0400 

0 

0 

0 

0 

0 

20 

20, 

1 

0 

0 

0 

0 

1 

o' 

0 

1 

0 

0 

0 

0 

C  = 

0 

0 

1 

0 

0 

D: 

= 

0 

0 

0 

0 

1 

0 

0 

0.001660 

0.514429 

0.000472 

0.001452  -0.068268 

0 

B2  F-16  QDES  Design  Information 

This  section  contains  the  feedback  gain  matrices,  Kalman  filter  gain  matrices,  and 
nominal  system  discrete  poles  for  the  various  (^ES  design  runs  for  the  F-16. 


B.2.1  F-16  1  g  Command 


K  =  [-0.01648  0.58854 

0.11327 

0.24277 

0.67653 

0.57277 

0  ■ 

0.71583 

0.50054 

0 

L  = 

0.72843 

0.47892 

0 

0.78094 

0.39017 

0 

0 

0 

0 

0.21629  ] 


Discrete  Poles: 
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0.41699±0.42016i 

0.28605 

0.74142 

0.88567 

0.98139 

0.99972  ±0. 003731 

0.99192 

0.28650 


B.2.2  F-I6  9.8  g  Command 


K  =  [-0.01648 


L  = 


0.58854 

0.11327 

0.24277 

0.97855 

0.08234 

0  ■ 

0.97913 

0.07223 

0 

0.98051 

0.06887 

0 

0.98735 

0.05445 

0 

0 

0 

0 

0.21629  ] 


Discrete  Poles: 


0.44808  ±0.541101 

0.28605 

0.74142 

0.88567 

0.98139 

0.99980±  0.00377i 

0.99192 

0.28650 


B.2.3  F-16 10.0  g  Command 

K  =  [-0.01648  0.58854  0.11327  0.24277  0.21629] 
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L 


Discrete  Poles: 


0.97873 

0.08071 

0 

0.97929 

0.07079 

0 

0.98066 

0.06750 

0 

0.98746 

0.05337 

0 

0 

0 

0 

0.4481010.541 16i 

0.28605 

0.74142 

0.88567 

0.98139 

0.99980±0.00377i 

0.99192 

0.28650 


B.2.4  F-16  10.2  g  Command 


K  =  [-0.01648  0.58854 

0.11327 

0.24277 

0.97890 

0.07913 

0  ■ 

0.97943 

0.06942 

0 

L  = 

0.98079 

0.06619 

0 

0.98757 

0.05233 

0 

0 

0 

0 

0.21629  ] 


Discrete  Poles: 


0.44812±0.54121i 

0.28605 

0.74142 

0.88567 

0.98139 

0.99980±0.00377i 

0.99192 

0.28650 


B.2J  F-16  10.4  g  Command 
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K  =  [-0.01648  0.58854 

0.11327 

0.24277 

0.97906 

0.07762 

0  ■ 

0.97957 

0.06809 

0 

L  = 

0.98093 

0.06493 

0 

0.98768 

0.05133 

0 

0 

0 

0 

0.21629] 


Discrete  Poles: 


0.44814±0.54127i 

0.28605 

0.74142 

0.88567 

0.98139 

0.99980±0.00377i 

0.99192 

0.28650 


B.3  F-16  Short  Period  Approximation 

The  state,  input,  and  output  vectors  with  dieir  associated  quantities  and  units  are  the 
same  as  for  the  A-4D  as  listed  in  Section  (A.3).  The  matrices  associated  with  the  model 
are  as  follows: 


■-1. 49100 

0.99600 

0.18800  ■ 

O' 

A  = 

9.75300 

-0.9600 

19.040 

B  = 

0 

0 

0 

0 

20 

1 

0 

0 

O' 

C  = 

0 

1 

0 

D= 

0 

_  0.54143 

0.00145 

-0.06827_ 

0 

B.4  F-16  Short  Period  Approximation  QDES  Design  Irtformation 


This  section  cmtains  the  feedback  gain  matrices,  Kalman  filter  gain  matrices,  and 
nominal  system  discrete  poles  for  the  various  QDES  design  runs  on  the  F-16  Short  Period 
approximation. 

B.4.1  F-16  Short  Period  I  g  Command 

K  =  [0.59479  0.21936  0.19813] 

‘0.71714  0.49985  0 

L=  0.78214  0.38947  0 

0  0  0  _ 

0.41735±0.42185i 

0.88681 

0.28605 

0.74049 

0.28650 

B.4.2  F-16  Short  Period  10.4  g  Command 

K  =  [0.59479  0.21936  0.19813] 

0.97958  0.06786  0 

L=  0.98770  0.05113  0 

0  0  0  _ 

0. 44840  ±0.54225i 
0.88681 
0.28605 
0.74049 
0.28650 

B.4.3  F-16  Short  Period  10.6  g  Command 

K  =  [0.59479  0.21936  0.19813] 
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0.97971  0.06659  0  ‘ 

L=  0.98780  0.05017  0 

0  0  0 

0.44841  ±0.542291 

0.88681 

0.28605 

0.74049 

0.28650 


Appendix  C.  MATLAB  Script  Files 

This  appendix  contains  representative  copies  of  the  MATLAB  script  files  written  by 
the  author  and  used  for  the  analysis  in  this  thesis.  They  are  recotnmended  as  models  for 
use  along  with  the  procedures  outlined  in  Chapters  V  and  VII. 


C.I  The ’'nzprp.m‘' File 

%  Filename:  "Qdes  souroefile  name  +  nzpip.m"  (for  noxnal  acceleration  prep) 

%  This  script  file  takes  the  aircraft  plant  data  for  the  longitudinal  axis 
%  and  fonns  the  lirpulse  Response  Matrices  (_tl,  _t2,  Jt3)  needed  by  Qdes  m 
%  well  as  the  jttt.mat  matrix  needed  to  form  the  controller.  To  use 

%  this  scr^,  you  need  the  cdrcraft  longitudinal  state  ^ce  plant  matrices 
%  and  the  actuator  c^namics.  You  will  also  need  to  have  the  nine  system 
%  matrices  deri^^  from  the  standard  P-K  form  of  the  aircraft  and  contoller 
%  set  up  for  NORMAL  AOCE1£RATION  (NZ)  command  %iith  your  desired  Exongenous 
%  Inputs  (w) ,  Regulated  Outputs  (z) ,  Sensed  Outputs  (y) ,  and  Actuator  Inputs 
%  (u) . 

%  Conments  in  ALL  CAPS  indicate  vtere  you  are  to  enter  data. 

%  Enter  plant  in  rad/  normalized  Velocity,  and  ft/sec^'Z  (Nz  at  eg) 

a= . 01266, - . 00588, - . 05075, 0; 

-. 101, -.8167,0,. 9984; 

0, 0, 0, 1; 

- . 21798, -12 . 4304, 0,-1.4219]; 

b=[0; 

-.08963; 

0; 

-19.388]; 


o=  (1,0, 0,0; 

0, 1, 0, 0; 

0, 0, 1, 0; 

0, 0, 0, 1; 

-64 . 034,  -517.788, 0,  -1 . 02708]  ; 

d=[0; 

0; 

0; 

0; 

-56.8254]; 

%  The  States  Are:  ubar  (normalized  forward  velocity) 
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% 

% 

% 

% 


2L4)ha  (angle  of  attack) 
theta  (pitch  angle) 
q  (pitdi  rate) 


% 

% 

% 

% 

% 


Outputs  are:  ubar 
alpha 
theta 

q 

nzog  (normal  g  at  the  center  of  gravity) 


%  Insert  actuator  dynamics 

ap*  [a,b;  zeros  (1,4),  -10] ; 

bp=  (zeros (4, 1) ; 10] ; 

cp“  [eye  (4) ,  zeros  (4, 1)  ;c (5, ; )  ,d(5, ; )  ]  ; 

dp*zeros  (5, 1) ; 


% 

% 

% 

% 

% 

% 


The  States  are  now:  ubar 
alpha 
theta 

q 

deltae 


NOTE:  All  subsequent 

<alcualtions  in  this  file 
as  well  as  the  nzpst.m  file 
that  prodLKses  the  results 
assime  these  states  in  this  order, 
(elevator  deflection) 


%  Outputs  are:  ubar 
%  alpha 

%  theta 

%  q 

%  nzog 

%  Change  sign  of  the  last  row  of  cp  (nzog)  to  make  positive  'g'  ip. 
cp(5,  :)=-cp(5,  :); 

%  Convert  to  deg,  ftsSec,  and  g's. 

T=>diag([634,57.3,57.3,57.3,57.3]);  %  State  Scaling 

S=[57.3];  %  Input  Scaling 

W=diag( [634,57.3,57.3,57.3, (1/32.2)]);  %  Output  Scaling 


^xl=T’*ap*inv(T) ; 
fcpcM'*bp*inv(S) ; 
qpd«W*qp*inv(T)  ; 
dp(>W*c^*inv(S) ; 

%  ENTER  THE  DESIRED  SCALING  FACR^.  This  allows  you  to  analyse  the  system 
%  for  other  than  a  unit  step  input.  If  you  put  in  put  in  1,  the  normal 
%  acceleration  ccnmand  iiput  is  1  g.  If  you  put  in  10,  the  normal  acceleration 
%  command  is  10  g's,  etc.  Note  that  all  the  constraints  are  still  writtei 
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%  as  if  you  axe  using  a  unit  step,  but  the  ou^xit  of  Nz  will  be  the  ze^xnse 
%  to  a  10  g  ocnnand  and  all  the  states  will  reflect  the  desired  10  g  oonnand. 
%  Note  that  you  MOST  also  enter  the  sann  scaling  factor  in  the  "nqsst.n^  file 
%  vtiich  produces  the  plots  of  the  system  response. 

Sfinz-1; 

cpds-  (l/Sfiiz)  *cpd; 

%  ENTER  THE  PLANT  DISCRETIZATIGN  RATE  IN  HZ. 


Ts-l/Srate; 

[aid,bd,ad,ckl]<>c2dm{apdi,bpd,cpds,dpd,Ts,  'zoh') ; 

%  EMIER  THE  NINE  dXS  STflC  MATRICES.  Note:  The  matrices  nust  be  labeled  (pa, 
%  pan,  pbu,  pcz,  pcy,  pdzw,  pdyw,  pdzu,  pdyu)  for  use  in  the  script  file 
%  "Qbc2ttt". 

pa«ad; 

pbi^[zeros(4,2);0, 10];  %  input  matrix  for  exogenous  inputs  (w) 

pbup*txl;  %  input  matrix  for  actioator  inputs  (u) 

pcz- [0,1,0, 0,0; 
od(5, :); 

0,0, 0,1,0; 

0,0, 0,0,1; 

zeros(l,5)];  %  output  matrix  for  regulated  outputs  (z) 

pcy- [0,0, 0,1,0  ; 
od(5, :) ; 

zeros  (1,5)];  %  ou^xjt  matrix  for  sensed  outputs  (y) 

pdz«^  [zeros  (4,2) ; 

0,1];  %  d  matrix  from  w  to  z; 

pdzu- [zeros (4, 1) ;1];  %  d  matrix  from  u  to  z 

pd^[0,0; 

0,0; 

1,0];  %  d  matrix  from  w  to  y 

pdyu-[0;0;0] ;  %  d  matrix  from  u  to  y;  zero;  plant  is  s.p. 

%  Form  an  observer  based  nominal  ccxitroller. 

%  Since  pydi^,  the  plant  is  strictly  proper.  Can  use  "obc2ttt".  Therefore 
%  need  "k"  the  State  Feedback  Matrix  and  "1"  the  Observer  Feedback  Matrix. 

%  Find  "k"  using  "dlqr" 

%  State  weighting  based  on  formula  (1/Xi''2)  where  Xi  is  max  value  the  state 
%  is  allowed  to  ta]ce  on. 

Ql^diag ( [ . 000044, . 0016, . 00111, . 00111, . 0016] ) ;  %  State  weighting 
Rl-.l;  %  Control  weighting 


C-3 


[k,  S,  E]  ^dlqr  (ad,  bd,  Ql,  Rl)  ; 


%  Find  "1"  Using  "dlqe" 

%  Put  "1"  on  every  state  that  oould  be  affected  by  a  plant 
%  disturbance  (wind  gust)  and  zero  on  others  (actuator  state) . 

Q2-.1; 

R2- .  001*qe  (size  (pd^,  1) )  ; 

/ 

l-dlqe(ad,G,pcy,Q2,R2)  ; 

%  Form  the  closed  loop  A  matrix  %fith  the  nominal  controller  ana  check 
%  its  poles. 

ACL>«[pa,-pbu*k;l*pcy,pa-l*pcy-pbu*k]  ; 

%  Use  ”obc2ttt"  to  get  the  Tl,  T2,  and  T3  Inpulse  Response  Matrices 

%  ENTER  THE  NAME  OF  THE  SOURCE  FII£  HERE 
%  I 

source  *  'A4D7NZ'; 

%  ENTER  THE  MICER  OF  SMf>l£S  This  is  the  plant  discretization  rate  (Srate) 

%  from  above  time  the  nunter  of  seconds  you  wish  to  analyse  the  system. 

%  (ex:  20  hz  times  10  seconds  200.  njsanple  ■*  200) 

n_sanple  »40  ; 

obc2ttt; 

dispC  '); 

disp ( ' ********************************************************************** ' ) ; 

dispC  '); 

di^CMien  writing  your  Qdes  source  file,  be  siue  that'); 
dispC  '); 

Zl= '  n_sattple“ ' ; 

Z2=  [Z1  nunCstr  (njsanple)  ] ; 

disp(Z2); 

dispC  '); 

disp  ('so  that  it  natches  with  the  inpulse  response  matrices  you  just  created'); 
dispC  '); 

Z8='  Hz  '; 

Z9='  sec'; 

Z5“'For  this  set  vp  your  sanple  rate  was:  '; 

Z7='Sanple  tine  was:  '; 

Z6-[Z5  nmiZstr  (Srate)  Z8  Z7  nizn2str(Ts)  Z9]; 

disp(Z6); 

dispC  '); 

dispC  '); 

dispCAlso,  be  sure  that  your  "nzpst.m"  file  uses  the  Scale  Factor'); 
dispC  '); 
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Z>'  Sfnz-  •; 

Z4-[Z3  nuRlZstr  (Sfhz)  ] ; 
clisp(Z4)  ; 
dispC  '); 
dispC  '); 

dispC  ********************************************************************** ' ) ; 
end; 


/ 

C.2  The  QDES  Source  File 

/*  Source  File  For  A4D7NZ  Filename:  "A4D'7NZ"  */ 

/*  This  Is  The  PREfiMBI£  */ 

/*  exogenous  inputs  (w)  */ 

♦define  NZJDCM  1  /*  reference  input  (Nz  ocmnand)*/ 

♦define  U_IN  2  /*  loop  irput  for  u  */ 

/*  regulated  outputs  (z)  */ 

1  /*  angle  of  attack  (alpha)  */ 

2  /*  nornal  acoel  (Nz)  */ 

3  /*  pitch  rate  (q)  */ 

4  /*  elevator  deflection  (deltae)  */ 

5  /*  elevator  deflection  ccmmanded  (deltaeccxm)  */ 

/*  sensed  outputs  (y)  */ 

♦define  Q_OUT_fffiAS  1  /*  measured  q  for  feedback  */ 

♦define  NZ_OUT_MEAS  2  /*  measured  Nz  for  feedback  */ 

♦define  NZ_BEF  3  /*  reference  output  (Nz  conmand)  */ 


/*  This  Is  The  DEJCLAFAIIONS  Section  */ 

/*  Note:  n_act*n_sens*n_tap  nwst  be  <  200  */ 

declarations  { 

njexog  =  2; 
n_reg  =  5; 
njsens  =  3; 
n_act  =  1; 
n_sanple  =  40; 
njtap  =  35; 

tl_coeffs  =  { 

♦include  "A4D7NZ_tl" 

}; 


♦define  AIFJXJT 
♦ctefine  NZjOUT 
♦define  QJDUT 
♦define  ELJDOT 
♦define  U  OUT 
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t2_ooeffs  ■>'  { 
♦inclucte  "A4D7NZ_t2« 

); 


tSjcoeffs  ■  { 
♦include  "A^7NZ_t3" 

}; 

) 


/*  Hiis  is  the  OCNSTOAINI^  Section  */ 
constraints  { 

/*  Put  COTStraints  On  U_CIUTA|_IN  Vector  Margins  */ 

for  w“0.0  to  %pi  st^  0.02  do  { 
n0g_H[5]  [2]  (1, w) <-1,6628; 

) 

/*  Put  Constraints  On  Elevator  Deflecticai  */ 

for  t=0  to  n_saitple-l  do 
-25.0<=u_st€p[4] (1]<=25.0; 

/*  Put  Constraints  on  Angle  of  Attack  */ 

for  t=0  to  n_sanple-l  do 
-5.0<=u_step[l] [1]<-18.0; 

/*  Put  Constraints  On  Overshoot  */ 

overshoot [2] [1]<=  0,1; 

/*  Put  Constrains  On  Undershoot  */ 

undershoot [2] [1]<=0.05; 

/*  Put  Constraints  On  Settling  Tine  */ 

for  t=16  to  20  do 

l-0.8392'^<=u_st^[2]  [1]  {t)<-l+0.8392'n;; 

/*  Put  Constraints  On  Steacfy  State  Error  */ 

for  t=20  to  njsanple-l  do 
0.97+0* (t-20)<=u_step[2] [1] (t)<=1.03-0*(t-20) ; 

/*  Put  Constraint  On  Rise  Time  */ 

for  t=12  to  n_sanple-l  do 
0.8<=u_st^[2]  [1]  (t); 

) 


C-6 


C3  The"mpstn"Fik 


%  Filenane:  "Qdes  source  file  name  +  nzpst.m”  (for  notnal  acceleration  post) 

%  Ihis  file  is  run  after  sucessfully  running  the  Qdes  source  file  and  getting 
%  the  Q  coefficients,  (files  .out  and  .q) 

%  This  file  takes  the  Qdes  fonned  Q  coefficients,  foms  the  Qdes  ccxitroller, 

%  forms  the  closed  loop  system,  and  siirulates  the  systan  respcme  to  a 
%  stqp  input  of  normal  acceleration.  It  outputs  a  plot  of  the  normal 
%  acceleraton  time  response  with  its  constraint  envelop,  angle  of  attack  and 
%  elevator  deflection  time  histories,  the  input  sensitvity  bode  magnitude  plot, 

%  the  closed  loop  bode  magnitude  plot,  and  the  input  vector  gain  and  plose 
%  margins. 

%  Note  that  if  you  used  a  scale  factor  in  the  "nzprp.m"  file,  the  plot  of  nontHl 
%  acceleration  will  be  a  scaled  to  represeit  the  response  to  your  desired 
%  commanded  normal  aooeleratic»i.  Ihe  time  histories  will  show  the  actual 
%  values  of  the  states  (angle  of  attack  and  elevator  deflecti<»i) . 

%  You  will  need  the  irputs  you  made  to  "visnz.m"  so  that  this  file  can 
%  draw  the  constraint  envelope  you  inposed  on  the  system  on  the  output  plot 
%  of  normal  acceleration. 

%  Remarks  in  ALL  CAPS  indicate  vdiere  you  are  to  enter  data. 

%  Use  "q6bc2oa(rp"  to  get  Q  coefficients  from  ".q"  and  form  the  controller 
%  Note:  "qpbc2ccrp"  uses  "jttt.mat"  Generated  in  "nzprp.m"  with  "obc2ttt" 

%  Therefore  it  nust  be  available  in  the  directory  you  are  worJdng  in. 

%  ENTER  THE  NAME  (X  THE  QXS  SOURCE  FII£  HERE 
%  I 

[ka,kb,kc,kd]=qc±)c2ccnp('A4D7NZ')  ; 

%  Use  "gethzw"  to  form  the  closed  loop  system  (Hzw) 

%  ENTER  THE  NAME  X  THE  QDES  SOURCE  FII£  HERE 

%  I 

[ha,  hb,  he,  hd]  =gethzw  ( '  A4D7NZ ' ,  ka,  kb,  kc,  kd)  ; 

%  Simulate  the  closed  locp  system  response  to  a  unit  step  input 

%  ENTER  PLANT  DISCRETIZATION  RATE  (IN  HZ)  AM)  TIME  YOU  WANT  THE  SMJIATICN 
%  TO  PUN  FCR.  Tliese  MUST  match  the  values  you  have  )3een  using  all  along  in 
%  "visnz.m",  "nzprp.m",  and  the  Qdes  source  file. 

Srate=8; 

Tmax=5; 

t=0:  (1/Srate)  :‘Hmax; 
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%  ENTER  THE  SCAI£  FACTOl  lOTl  IN  "nzprp.m" 

Sfinz=l 

hC3-[hc(l, :)  ;Sfin2*hc(2, :)  ;hc(3, :)  ;hc(4, :)  ;hc(5,  :)]; 

%  Create  the  input  vector  (w  for  the  closed  loop  system  Hzw) 

%  ENTER  INPUTS  EXB  THE  SIMULATION.  NOTE  THAT  THE  CBIKR  OF  TOE 
%  INPUTS  MUST  MATCH  THAT  USED  IN  THE  PREAEBIE  OF  THE  CPES  SOURCE  FII£. 

ul*ones (length (t),l);  %  Refereioe  (commanded  normal  acceleration)  (always  1) 

u2=zeros (length (t) ,1) ;  %  U_IN  (always  zero) 

u=[ul,u2]; 

[y,X]=dlsim(ha,hb,hcs,hd,u) ; 

%  ENTER  DATA  FROM  "visnz.m"  FOR  THE  NC«MAL  AOCEIERATION  CONSTRAINT  ENVELOPE 

Tmz=1.5; 

Hxiz=0.1; 

Usnz=-.05; 

Tsnz=2.5; 

pemz^S; 

Es3n2=3; 

Tsnzst=2.0; 

Ynz=.2; 

Thpnz=1.5; 

Tyznz=0.5; 

Tynznz=1.5; 

ytynznz=.80; 

%%%%%%%%%%%%  calculations  for  nontal  acoeleratiOT  envelope  plot  %%%%%%%%%%%%%% 

tlnz=0;t2nz=0;t3nz=0;t4nz=0;t5nz=0;t6nz=0;t7nz=0;t8nz-0;npnz=0;tmz=0; 
t  sn2up=0;  tsnzlo=0 ;  essnzi;p=0 ;  essnzlo=0 ;  hil  30^=0 ;  hi2  3nz=0 ;  lol3nz=0 ;  lo2  302=0 ; 
usnz=0; 

Stime=l/Srate; 

E^z=Mpnz+l; 
t lnz=0 : St ime : Tsnzst ; 
for  i=l :  length  (tlnz) , 
itpnz  (i)^^z; 

end; 

npnz=^ipnz*Sfnz  ; 

for  i=l :  length  (tlnz) , 
usnz(i)^snz; 

end; 

usnz=usnz*Sfnz; 


%  Pea]c  Overshoot 


%  Ifedershoot 
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%  Rise  Time 


tZnz^Tmz ;  St  iire ;  'Dnax; 
for  i=l : length (t2nz) , 
tmz(i)=0.8; 

end; 

tmz=tmz*Sfnz; 

t3nz=^3nzst:StinE;Tsnz;  %  Settling  Time 

Tsnzeb=( .  01*pemz)  ^  (1/  (Tsnz*Srate) )  ; 
for  i=l :  length  (t3nz) , 

t3nzup(il“l+Tsnz^''  <i+Tsnzst*Srate)  ; 
tsnzlo  (i)  =1-T3nzeb'‘  (i+T3nzst*Srate)  ; 

end; 

tsnzi;p=tsnzi;p*Sfiiz  ; 
tsnzlo=t3nzlo*Sfinz ; 

if  Essnz==0,  %  Steac^  State  Error  When  Set  To  Zero 

t4nz=0:Stime:lhBx;  %  Cc»itinues  Curve  Fran  Settling  Time 

for  i=T3nz*Srate:  ( ('nTiax*Srate)  +1) , 
esanzijp  (i  =1+Tsnz^''  (i)  ; 
essnzlo  (i)  “l-Tanz^"  (i) ; 

end; 

else 

t4nz=Tsn2:Stime:'Ihax;  %  Steacfy  State  Error  When  Not  Zero 

%  Draws  Line  Fran  Y  Value  0  Tsq  To  Essq  0  Unax 

essnzdcr=  ( ( 14pemz* .  01)  -  ( 1+Essnz* .  01) )  /  (!Ilnax*Srate-T3nz*Srate) ; 
for  i=l :  ( (linax*Srate-T3nz*Srate)  +1) , 

esanzup  (i)  =  (1-fpemz*  .01)  -essnzdcr*  (i) ; 
essnzlo  (i)  =  (1-pemz* .  01)  4essnzdcr*  (i) ; 

end; 

end; 

essnzvp=essnzi5J*Sf nz ; 
essnzlo=essnzlo*Sfnz ; 

hillnz=(  .5* (J^pnz-Ynz) ) ''  (1/3)  ; 
hil2nz=hilln2/ ( . 4*Dtpnz*Srate) ; 
t5nz=0 :Stime :  ( .  4*Titpnz)  ; 
for  i=l:  (.4*'IVtpnz*Srate)+l, 

hil3nz  (i)  -  (hil2nz*  (i-1) )  ''3+Ynz; 

end; 

hil3n2=hil3nz*Sfnz; 

hi21nz=  ( (^%3nz- .  5*  (H^nz-Ynz) )  -Ynz) "  (1/3) ; 
tdi£hinz=  (Thpnz*Srate)  -  ( .  4*Thpnz*Srate) ; 
hi22nz=hi21nz/tdifhinz; 

1 6nz= .  4  *Thpnz :  St  ime :  Tlrrpnz  ; 
for  i=l:tdifhinz+l, 

hi23n2  (i)  =  ( ( (^^z- .  5*  (t^z-Ynz) )  -Ynz)  -  ( (hi21nz-hi22nz*  (i-1) )  ''3) )  ... 

+  (^^z- .  5*  (f^z-Ynz) ) ; 

end; 

hi23nz=hi23nz*Sfnz; 
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lollnz”  ( .  5*Ytynzn2)  ( 1/3) ; 

tdif  lolnz= .  4*  { (Tynznz-Tyznz)  *Srate)  ; 

lol2n2=lollnz/tdiflolnz; 

t7nz=Tyznz: Stine: (Tyznz+.4* (Tynznz-Tyznz) ) ; 

for  i=l:tdiflolnz+l, 

lolSnz  (i) » (lol2nz*  (i-1) )  ''3; 

end; 

lol 3nz=lol3nz *Sf  nz ; 

lo21nz=(  .5*  ^ynznz) ) ''  (1/3) ; 
tciiflo2nz=. 6*  ( (Tynznz-Tyznz)  *Srate) ; 
lo22nz=lo2 Inz/tdif lo2nz ; 
t8nz=. 4* (Tynznz-Tyznz) +Tyznz : Stine  rTynznz; 
for  i=l:tciiflo2nz+l, 

lo23nz  (i)  = .  5*  (Ytynznz)  -  (lo21nz-lo22nz*  (i-1) )  '^3+ .  5*Ytynznz; 

end; 

lo2  3nz=lo2  3nz*Sfnz ; 


%%%%%%%%%%%%%  plot  nomel  g  tine  response  and  omistraint  envelope  %%%%%%%%% 
figured) ; 


subplot (111) ; 
plot(t,Y(:,2)); 

%title('A4D  NZ_OUT/NZ_OCM'); 
ylabel  ( 'Nomal  g  (gs) ' )  ; 
xlabeK'Tine  (sec)'); 
grid; 
hold; 

plot  (tlnz,npnz, ' : ' )  ; 
plot  (tlnz,  usnz, ' : ' )  ; 
plot  (t2nz,tmz, ' : ' ) ; 
plot (tSnz/tsnztp, ' : ' ) ; 
plot (t3nz,tsnzlo, ' : ' ) ; 
plot  (t4nz,es3nzi4), ' : ' )  ; 
plot (t4nz,e3snzlo, ' : ' ) ; 

%plot  (t5nz,hil3nz, ' : ' )  ; 

%plot (t6nz,hi23nz, ' : ' ) ; 

%plot (tlnz, lol3nz, ' : ' ) ; 

%plot  (tSnz,  lo23nz, ' : ' )  ; 
hold; 


%  plots  normal  g  tine  response 


%  plots  peak  overshoot  line 
%  plots  undershoot  line 
%  plots  rise  tine  line 
%  plots  iQiper  settling  tine  curve 
%  plots  lower  settling  tine  curve 
%  plots  v?:per  steady  state  error  curve 
%  plots  lower  steady  state  errorcurve 
%  plots  1st  half  of  i??«r  constraint  curve 
%  plots  2nd  half  of  upper  constraint  curve 
%  plots  1st  half  of  lower  constraint  curve 
%  plots  2nd  half  of  lower  constraint  curve 


t  lnz=0 ;  t2nz=0 ;  t3nz=0 ;  1 4nz=0 ;  1 5nz=0 ;  1 6nz=0 ;  t7nz=0 ;  1 8nz=0 ;  npnz=0 ;  t  mz=0 ; 

tsnzup)=0;tsnzlo=0;essnztp=0;essnzlo«0;hil3nz=0;hi23nzK);lol3nzK);lo23nz=0; 

usnz=0; 

%%%%%%%%%%%%%%  plot  angle  of  attacdc  and  elevator  deflection,  %%%%%%%%%%%%% 

figure (2) ; 

sutplot  (211) ; 
plot(t,Y(:,4)) ; 
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%title( 'Elevator  Deflecticxi  to  Noimal  G  Ccnmand*)  ; 

ylabel ( 'deltae  (deg) ' ) ; 

grid; 

sufcplot (212) ; 
plot(t,Y(:,l)) 

%title( 'Angle  of  Attack  to  Normal  G  Caatard'): 
ylabel ( ' alpha  (deg) ' ) ; 
xlabelCTime  (sec)'); 
grid; 

%%%  figure  M-C^ircle  radius  and  gain  aivl  phase  margins  for  U_CUrAJ_IN  %%%%% 

[M,  P,  W]  =cfcode  (ha,  hb,  he,  hd,  Stime,  2) ; 
m_crcl=l/ (max (M ( : , 5) ) ) ; 
upgm(±i^-20*logl0  (l-m_crcl) 
lwgtndti=-20*logl0  (l-rtn_crcl) 
ptrde^57 . 3*2*asin  (m_crcl/2) 

%%%%%%%%%%%%%%%%%%  plot  ir^xit  sensitivity  of  U_OUTAJ_IN  %%%%%%%%%%%%%%%%%% 

M®=20*logl0  (M)  ; 

figure (3) ; 

subplot (211) ; 

semilogx (W,M( : ,  5) )  ; 

ylabel ( 'Magnitude  (not  dB) ' ) ; 

grid; 

%%%%%%%%  plot  conplimentary  sensitivity  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

tdb^O; 

[M2,  P2,  W2  ]  =dbode  (ha,  hb,  he,  hd,  Stime,  1) ; 

M2dB=20*logl0(M2); 

sufcplot (212) ; 

semilogx  (M2,M2dB  ( : ,  2) ) ; 

grid; 

hold; 

for  i=l : length (W2) , 
tdb(i)=-3; 

end; 

semilogx  (VC,  tdb, ' : ' ) ; 

%title( 'Conplimentary  Sensitivty  NZ_OOT/NZ_OCM' ) ; 

ylatel  ( 'Magnitude  (dB) ' )  ; 

hold; 

xlabel ( 'Frequency  (rad/sec) ' ) ; 

%%%%%%%%%%%%  output  the  size  of  the  Qdss  constroller  %%%%%%%%%%%%%%%%%%%%%%%% 
orc^size  (ka,  1) ; 


C-11 


Xl-^um23tr  (oid)  ; 

X2“’The  Order  Of  The  Qdes  Controller  is 
Z1-IX2  XI]; 
dispC  '); 
dispC  '); 

Hi  !Rp  ^  *  **AA**A*ikikA***1hA*ik*A***A**i^*ik*A********4k**1kilrik**ikrik******A**ik****  f  ^  ■ 

dispC  '); 
disp(Zl)  ; 
dispC  '); 

dli>?p(  '  ***^**A*A*** **************** *******************************11:***  >  ^  • 
end; 


C.4  The"vi  zm"  File 

%  This  file  assists  in  developing  the  unit  step  response  cx>nstraint  envelope 
%  for  the  aircraft  noxnal  acceleration  (Nz) 

%  Filename;  "visnz.m". 

%  In  addition  to  dravdng  the  containt  envelope,  the  program  generates  the 
%  Qdes  Source  File  statements  (in  proper  fomat  and  syntax)  to  force  Qdes  to 
%  iiipose  the  exact  same  constraint  envelope  as  that  developed  by  this  program. 

%  Also  the  program  will  generate  the  Qdes  Source  File  statements  that  iirpose 
%  vector  margins  of  your  choice  on  any  desired  SISO  transfer  function  of  the 
%  system. 

%  The  program  requires  that  you  enter  the  desired  Rise  Time  (Tr) ,  Settling  Time 
%  (Tsl)  and  desired  percent  (per) ,  Overshoot  ()^) ,  ttidershoot  (Us) ,  and 
%  Steady  State  Error  (Ess)  for  the  time  re^x5nse  of  Nz. 

%  Addtional  constraint  barriers  can  be  constnocted  that  do  not  correspond  to 
%  any  particular  stqp  response  figure  of  merit  )3ut  are  useful  in  shaping  t)ie 
%  rqponse.  These  are  described  prior  to  the  section  «here  you  are  required 
%  to  input  the  ^»cificatiOTS . 

%  The  program  also  requires  you  to  enter  the  discretization  parameters  that  you 
%  used  vhen  setting  vp  the  control  problem  for  Qdes.  This  involves  the 
%  Sanple  Rate  (in  Hz)  and  the  time  period  you  asJced  Qdes  to  wor)c  over  (in  sec) . 
%  This  turns  out  to  be  "n_sanple"  times  the  sanple  time  which  is  the  inverse 
%  of  the  Sanple  Rate  (ie.  1/20  Hz  =  .05  secs) . 

%  Carments  in  ALL  CAPS  indicate  \*iere  you  are  to  enter  data. 

%  ENTER  THE  DESIRED  VBCTC»  MARGINS  FOR  THE  SISO  TRANSFER  FUNCTION  TOU  WISH  TO 
%  CONSTRAIN.  LOWER  GAIN  MARGIN  MUST  BE  NEGATIVE. 

Dugnnz=7;  %  Desired  Upper  Gain  Margin  (dB)  Range:  (0.0  <  Dugnnz  <  +inf) 

Dlgrnnz=-4;  %  Desired  Lower  Gain  Margin  (dB)  Range:  (-6.0  <  Dlgnnz  <  0.0) 

Dpmnz=35;  %  Desired  Phase  Margin  (deg)  Range:  (0.0  <  I^mnz  <  60) 

%  ENTER  TIME  PERIOD  (IN  SECS)  PC»  WHICH  YOU  ASKED  CPES  TO  WC»K 
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%  ENTER  SMflE  RATE  (IN  HZ)  AT  WHICR  YOU  DICRITIZED  THE  PLANT 
Srate-8; 

Stiine=l/Srate; 


%  ENTER  THE  FIVE  UNIT  STEP  FIGURES  (F  MERIT  YOU  WISH  GPES  TO  MEET 


Tmz“1.5; 

Hxiz“0.1; 

Usnz=-0.05; 

Tsnz=2.5; 

pemz*3; 

Essnz=3; 


%  rise  tiite 
%  overshoot 

%  undershoot  (enter  negative  nurrber) 

%  settling  time 

%  percent  value  at  the  settling  time 
%  steacfy  state  error  (most  be  <=■  pemz) 


%  The  Settling  Time  constraint  envelope  is  an  ejiponential  curve  that  starts 
%  at  a  time  of  your  choosing,  decays  e^qpcai&itially  while  it  passes  throu^ 
%  the  proper  percentage  above  and  below  1  at  the  requested  Ts,  and  achieves 
%  the  requested  Ess  at  Trrax.  (This  is  why  Ess  must  be  <-  pemz) . 


%  ENTER  THE  TIME  YOU  WISH  THE  SETTLING  TD©  OONSTRAINr  TO  BECIN 


Tsnzst=2; 

%  You  may  refine  the  i^^per  constraint  envelc^je  with  an  e^qponetial  curve  >*iich 
%  starts  at  a  point  of  your  choosing  cn  the  y  axis  at  time  zero,  and  meets  the 
%  overshoot  line  at  a  time  of  your  choosing. 

%  ENTER  Y-VAIUE  AT  TIEE  ZERO  WERE  YOU  WANT  THE  UPPER  (Y^ISTRAINr  LINE  TO  START 
Ynz=.2; 

%  ENTER  THE  TIME  AT  WHICH  YOU  WANT  THE  UPPER  CONSTRAINT  LINE  TO  TOUCH  THE 
%  OVERSHOOT  LINE 


'nTpnz=1.5; 


%  You  may  refine  the  lower  constraint  envelcpe  with  an  e;gx5nential  curve  that 
%  starts  at  y=0  at  a  tiite  of  your  choosing  and  ends  at  a  Y  value  and  time 
%  of  your  choosing. 

%  ENTER  THE  TIE©  WHERE  YOU  WANT  THE  ICSttR  CCNSTRAIOT  ENVELWE  CURVE  TO  START 


'iyznz=0.5; 
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%  ENTER  THE  POINT  (TIME  AND  Y  VAUIE)  tiiERE  YOU  »iANT  THE  UMER  OCNSIRMNT 
%  ENVELOPE  CURVE  TO  END  (noimally  at  0.8  and  Tr) . 

Tynznz=l . 5; 

Ytynznz=.80; 

%%%%%%%%%%%%%%%%%%%%%%  calculaticxis  for  nozinal  g  %%%%%%%%%%%%%%%%%%%%%%%%%%% 

tln2=0;t2nz=0;t3nz=0;t4nz=0;t5nz=0;t6nz»0;t7n2=0;t8nz=0;npnz*0;tmz*0;tsimp'=0; 
tsnzlo=0;essnzup^0;essnzlo*0;  hil3nz=0;hi23nz=0;  Iol3n2=0;  lo23nz=0/usnz=0; 

Senslnz*!/  (1-10''  (-Dugptnn2/20) ) ; 

Sens2nz=l/  (10''  (-Dlgimz/20)  -1) ; 

SensSnz®!/  (2*sin  ( (E|xmz/57 .2958)  /2) ) ; 

Sens4nz°°[Senslnz,Sens2nz,Sens3nz] ; 

Sensnz^nin(Sens4nz)  ; 

^anz=4^z+l; 
t  lnz=0 :  St  lire :  Tsnzst ; 
for  i=l : length (tlnz) , 
npnz(i)^4xiz; 

end; 

for  i=l : length (tlnz) , 
u3n2(i)=U3nz; 

end; 

t2nz*Tmz :  Stime :  Thax; 
for  i=l: length (t2nz), 
tmz(i)=0.8; 

end; 


%  Peak  Overshoot 


%  Undershoot 


%  Rise  Time 


t3nz=T3nzst:Stiine:Tsnz;  %  Settling  Time 

Tsnzeb=  ( .  01*pemz) ''  (1/  (Tsnz*Srate) ) ; 
for  i=l :  length  (t3nz) , 

tsnzup  (i)  =1+Tsnz^''  (i+Tsr.z3t*Srate) ; 
t3nzlo  (i)  =1-Tsnzeb''  (i+Tsn23t*Srate) ; 

end; 


if  E33n2=0,  %  Steacfy  State  Error  When  Set  To  Zero 

t4nz=0:  Stime  :Tltax;  %  Continues  Curve  Fran  Settling  Time 

for  i==T3nz*Srate:  ( (T)iBX*Srate)  +1) , 
essnzif)  (i)  =1+T3nzeb''  (i) ; 
essnzlo  (i)  =1-T3nzeb''  (i) ; 

end; 

else 

t4nz=Tsnz:  Stime  rUnax;  %  Steady  State  Error  When  Not  Zero 

%  Draws  Line  Fran  Y  Value  0  Tsq  To  Essq  @  Dnax 

essnzdcr=  ( (I4pemz* .  01)  -  (1+Essnz* . 01) )  /  (Thiax*Srate-T3nz*Srate) ; 
for  i=l :  ( (TSnax*Srate-Tsnz*Srate)  +1) , 

essnzv?)  (i) = (14pemz* .  01)  -essnzcfcr*  (i) ; 
essnzlo(i)=(l-pemz*.01)4essnzdcr*(i) ; 
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end; 

end; 


hillnz-  ( .  5*  (Wlpnz-Ynz) ) ''  (1/3)  ; 
hil2nz“hillnz/  ( .  4*Tlqpnz*Srate) ; 
t5nz=<0:Stiine:  (.4*lhpnz) ; 
for  i»l :  ( .  4*'nr{anz*Srate)  +1, 

hil3nz  (i)  -  {hil2nz*  (i-1) )  ''3+Ynz; 

end; 


hi21nz“  ( (J^jnz- .  5*  (Hpnz-Yhz) ) -IDnz) '' (1/3) ; 
tdifhin^  (Tlt({»z*Srate)  -  ( .  4*Tlqpnz*Srate) ; 
hi22nz*4u21nz/tdl£hinz; 
t&iz> .  4*'Dipnz :  Stine  iTkipnz; 
for  i=l:tdifhinz+l, 

hi23nz  (i)  -  ( ( (tifnz- .  5*  (H»z-Yhz) )  -»«;)  -  ( (hi21nz-*i22nz*  (i-1) )  *3) )  ... 
+  (»t«z- .  5*  a^-Ynz) )  ; 

end; 


lollnz- ( . 5*Ytynznz) " (1/3) ; 

tdif lolnz” . 4* ( (Tynznz-Tyznz) *Srate) ; 

lol2nz=lol Inz/tdif lolnz ; 

t7nz=Tyznz: Stine: (Tyznz+,4* (Tynznz-Tyznz) ) ; 

for  i«l:tdiflolnz+l, 

lol3nz (i)-(lol2nz* (i-1) ) ^3; 

end; 

lo21nz=‘  ( .  5*  (Ytimznz) )  (1/3) ; 

tdiflo2nz*. 6* ( (Tynznz-Tyznz) *Srate) ; 

lo22nz=lo21nz/tdiflo2nz; 

t8nz= . 4* (Tynznz-Tyznz) +Tyznz : Stine  :Tynznz; 

for  i=l:tdiflo2nz+l, 

lo23nz (i) = . 5* (Ytynznz) - (lo21nz-lo22nz* (i-1) ) "3+ . 5*Ytynznz; 

end; 


%%%%%%%%%%%%%%%%%%%  plot  nomel  g  constraint  envelope  %%%%%%%%%%%%%%%%%%%%%%% 


figure (1) ; 

sulplot  (111) ; 

plot  (tlnz/npnr, ' : ')  ; 

%  plots 

hold; 

plot (tlnz,usnz, ' : ' ) ; 

%  plots 

plot  (t2nz,tmz, ' :  ’ ) ; 

%  plots 

plot (t3nz,tsnzvp, ' : ' ) ; 

%  plots 

plot (t3nz,tsnzlo, ' : ' ) ; 

%  plots 

plot  (t4nz,essnzup, ' : ') ; 

%  plots 

plot  (t4nz,essnzlo, ' : ' ) ; 

%  plots 

plot (t5nz,hil3nz, ' ; • ) ; 

%  plots 

plot (t6nz,hi23nz, ' : ' ) ; 

%  plots 

plot (t7nz, lol3nz, ' : ' ) ; 

%  plots 

plot (tSnz, lo23nz, ' : ' ) ; 

%  plots 

hold; 

grid; 

peaJc  overshoot  line 

undershoot  line 
rise  tine  line 
vpper  settling  tine  curve 
lower  settling  tine  curve 
vpper  steady  state  error  curve 
lower  steady  state  error  curve 
1st  lialf  of  upper  constraint  curve 
2nd  half  of  upper  ccnistraint  curve 
1st  half  of  lower  ccxistraint  curve 
2nd  half  of  lower  ccxistraint  curve 
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axis  { I0,ltaax,Usn2-.l,»%*Jz+.2] )  ; 
xlabelCtiroe  (secs)'); 
ylabel  <  'noxnal  g  (Nz) ' ) ; 

title  ( 'Noznel  G  (Nz)  Unit  St^  Re^xxise  (Constraint  Envelope ' ) ; 

tlnz«0;t2nz^;t3nz^;ttoH);t&iz^;t6nz^;t7nz^;t8nz=0;npnz“0;tmz«0;tsnzupN); 

t8nzlo“0;essnzvjqp^;essnzlcH);hil3nz-0;hi23nz-0;lol3nz=0;lo2^iz*0;usnz-0; 

%%%%%%  generate  Qcles  source  file  constraint  statements  %%%%%%%%%%%%%%%%%%%%%% 


Al-'oiwershootlout]  Iin]<“ 
A2“'for  t='; 


A3='  to 
A4='  do'; 

A5“*n_sanple-1  do'; 
A6“'<=u_step[out]  tin]  (t)<»'; 

A8='l+'; 

A^'l-'; 

A10='+'; 

All='*'; 

A12=' ('; 

A13=')  '; 


A15='t'; 

A16='-'; 

Sl=';'; 


Al^*u_ateptout]  [in]  (t)<='; 
A2(^'"3'; 


A25-'<“u_step[out]  [in]  (t) ;  '; 


A31“'for  »^.0  to  ^i  step  0.02  do  ['; 
A32“'mag_H[out]  [in]  (l,w)<-'; 

A3>')'; 

A34-'undershoot[out]  [in]<>*'; 


dispC  '); 
disp  ( '  ' ) ; 

dlsp  ( *  *********T******************'^**********^******************** *********** * )  ; 
( * ********************************************************************** * )  ; 
disp ( '  ' ) ; 

disp ('This  is  how  your  Qdes  Source  File  constraint  statements  should'); 
dispC  '); 

dispClocJt  to  draw  the  equivalent  constraint  envelope  for  Qdes  and  oraistrain'); 
disp ( '  ' ) ; 

disp  ('the  vector  margins  of  your  selected  SISO  transfer  function. '); 
dispC  ’); 


%%%%%%%%%%%%%%%%%%%%%%%  for  normal  acceleration  %%%%%%%%%%%%%%%%%%%%%%%%%%%% 


A14n=num2str  (essnzdcr) ;  A26n=nuii2str  (Ytynznz* .  5)  ; 

A17n=num23tr  (1-pemz*  .01) ;  A27n=TOnCstr  (lo21nz)  ; 

A18n=niJti2str  (14pem2* . 01) ;  A28n=num2str  (lo22nz) ; 

A21n=num2str  ( (Haaz- .  5*  (Hpnz-Ynz) )  -Ynz)  ; 

A22n^um2str  (^xiz-.  5*  (Hanz-Ynz) ) ; 

A23n=nura2str(hi21nz) ;  A29n=nun2str(Srate*(TVznz+.4*(TVnznz-TVzr-r) ) )  ; 

A24n=num2str  (hi22nz) ;  A30n-inan2str  ( .  4*'Ilrpnz*Srate) ; 

Tln=num2str  (Tmz*Srate) ; 
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T2nrmin2str  (Tsnz*Srate)  ; 

T3n-nwni2str  (Tanzst^Srate)  ; 

T4n'4iutCstr  C&ipnz^Srate) ; 

T&i-niin2str  (lyznz^Srate) ; 

TCnpnumZstr  (Tynznz*Srate) ; 

Bln-[A1  num2str(^|anz-l)  SI]; 

B2n-[A2  T3n  A3  T2n  A4]; 

B3n“[A9  nun23tr (Tsnz^)  A7  A6  Afl  nun2str (Tsnzeb)  A7  SI]; 

B4n-[A2  T2n  A3  A5]; 

B5n-[A17n  AlO  Alta  All  A12  A15  A16  T2n  A13  A6  A18n  A16  Alta  All  A12  A15  A16  ... 

T2n  A13  SI]; 

B6n»[A2  Tin  A3  A5]; 

B7n-[A2  nun2str(0)  A3  A30n  A4]; 

B8n-[A19  A12  nuiCatr (hil2nz)  All  A15  A13  A20  AlO  num2atr(Ynz)  SI]; 

B9n-[A2  A30n  A3  Tta  A4]; 

BlQn-[A19  A21n  A16  A12  A23n  A16  A2ta  All  A12  A15  A16  A3Qn  A13  A13  A20  AlO  ... 
A22n  SI]; 

Blln=[A2  T5n  A3  A29n  A4]; 

B12n=[A12  nvjn23tr(lol2nz)  All  A12  A15  A16  T&J  A13  A13  A20  A25]; 

B13n-[A2  A29n  A3  T6n  A4]; 

B14n>[A26n  A16  A12  A27n  A16  A28n  All  A12  A15  A16  A29n  A13  A13  A20  AlO  A26n  A25]; 
B15n*[A32  nunCstr (Sensnz)  SI]; 

Bl&v-[A34  nLBn23tr(-U3nz)  SI]; 

dLap ('*****»*********  EX»  tCTMRL  MCEIERATICN  <Nz)  **********•); 
diapC  *); 

diapC/*  Put  Conatrainta  On  "  "  Vector  Margina  */•); 

diapC  '); 

di3p(A31)  ; 

di^(B15n) ; 

di3p(A33) ; 

diapC  '); 

diapC/*  Put  Conatrainta  On  Overahoot  */'); 
di^C  '); 
diaptBln) ; 
di^C  '); 

diapC/*  Put  Ccxiatraina  On  UtxJerahoot  */'); 
diapC  '); 
di3p(B16n) ; 
diapC  '); 

di3p(’/*  Put  Conatrainta  On  Settling  Time  */'); 

diapC  '); 

di8p(B2n) ; 

disp(B3n) ; 

diapC  '); 

di^C/*  Requeat  Zero  Steady  State  Error*/'); 
di^C  '); 

di3pCRe_H[out]  [in]  (1,0)  =1;'); 
diapC  '); 
diapC  (»'); 
diap  C  ' ) ; 
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<3xsp{’/*  Put  Constraints  On  Steac^  State  Error  (if  not  zero)  */'): 

di^C  '); 

di^(B4n)  ; 

di^  (B5n)  ; 

di^C  '); 

diqpC/*  Put  Constraint  On  Rise  Time  */’); 
dispC  '); 
disp(B6n)  ; 

di2p(’0.8<=nj_step[out]  [in]  (t) ;  •); 
dispC  '); 

dispC/*  Put  In  Upper  Contraint  Curve  (if  desired)  */’); 

di^C  '); 

cLLsp(B7n) ; 

disp (B8n) ; 

dispC  '); 

disp (B9n) ; 

disp (BlOn) ; 

dispC  '); 

dispC/*  Put  In  liower  Constraint  Curve  (if  desired)  */'); 

dispC  •); 

disp(Blln); 

di^  (B12n)  ; 

dispC  '); 

di^  (B13n) ; 

disp(B14n) ; 

dispC  •); 

dispC  •); 

cid.3p(  *  **★★**•***•*!•***♦■*** **'fc*****-*'*l*li*'fc*'*'fc^*'A**'<Ht******^**'*'fc<i*'*'fc'<i-*'J<1»****  • 

dispC  •); 
dispC  •); 
flpq=*lpq-i; 

H^nz-t^anz-l; 

disp ('Here  are  the  nutbers  to  put  into  the  "nzpst.m"'); 

di^C files  to  draw  the  tine  sane  tine  response  curves  you  just  created'); 

dispC  •); 

dispC***************  FOR  NOWftL  AXEIfRATION  (Nz)  *****★****'); 
dispC  '); 

Tmz 

F(^z 

Usnz 

Tsnz 

pemz 

Essnz 

Tsnzst 

Ynz 

TSrpnz 

Tyznz 

Tynznz 

Ytynznz 

disp ( '**********************************************************************' ) ; 

J  • ★**★***★*★★*★****★*******★*★★*************★**♦**★★*★*****★**★****★★★** • J  • 
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end; 


CJ  The  "surf m"  File 

%  This  program  calcxilates  the  control  svirfaoe  deflection  rate  from  the  time 
%  history  of  the  control  siarfaoe 

%  To  use,  adjust  plot  axes  so  that  you  can  ID  the  steepest  slope  of  the 
%  control  surface  deflecticm  time  history.  Use  the  cursor  to  pull  off  the 
%  nagnitude  and  tine  of  the  two  end  points  of  the  segment  vdth  the  largest 
%  slope.  The  program  will  taJce  the  absolute  value  of  the  difference  of  the 
%  magnitixles  and  devide  it  by  the  tine  difference  between  the  two  to  get  the 
%  control  surface  deflection  rate. 

[tl,ml]=ginput  (1)  ; 

[t2,m2]=girput  (1) ; 

mdif  f=abs  (m2-ml)  ; 
tdif f=abs (t2-tl) ; 

max_deflection_rate=mdiff/tdiff 

%  Set  varibles  back  to  zero 

tl=0;  t2»0;  ml=0;  m2=0;  mdiff=«0;  tdiff=0;  max  deflection  rateK); 


end; 
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